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ABSTRACT 

Deep optical observations of the spectra of 12 Galactic planetary nebulae (PNe) and 3 Mag- 
^> ' ellanic Cloud PNe were presented in Paper I by Tsamis et al. (2003b), who carried out an 

! abundance analysis using the collisionally excited forbidden lines. Here, the relative inten- 

OO . sities of faint optical recombination lines (ORLs) from ions of carbon, nitrogen and oxygen 

are analysed in order to derive the abundances of these ions relative to hydrogen. The relative 
intensities of four high-^ C II recombination lines with respect to the well-known 3d-4f A4267 
line are found to be in excellent agreement with the predictions of recombination theory, re- 
moving uncertainties about whether the high C'^^ abundances derived from the A4267 line 
could be due to non-recombination enhancements of its intensity. 
i-C ' We define an abundance discrepancy factor (ADF) as the ratio of the abundance derived 

for a heavy element ion from its recombination lines to that derived for the same ion from its 
Q ■ ultraviolet, optical or infrared collisionally excited lines (CELs). All of the PNe in our sample 

}-i ■ are found to have ADF's that exceed unity. Two of the PNe, NGC 2022 and LMC N66, have 

, ADF's of 16 and 11, respectively, while the remaining 13 PNe have a mean O^^ ADF of 

■ 2.6, with the smallest value being 1.8. 

Garnett & Dinerstein (2001a) found that for a sample of about a dozen PNe the magnitude 
of the 0^+ ADF was inversely correlated with the nebular Balmer line surface brightness. We 
have investigated this for a larger sample of 20 PNe, finding weak correlations with decreasing 
I surface brightness for the ADF's of 0'^+ and C"^^. The C'^^ ADF's are well correlated with the 

absolute radii of the nebulae, though no correlation is present for the ADF's. We also find 
both the and ADF's to be strongly correlated with the magnitude of the difference 
between the nebular [O III] and Balmer jump electron temperatures (AT), corroborating a 
result of Liu et al. (2001b) for the ADF. AT is found to be weakly correlated with 
decreasing nebular surface brightness and increasing absolute nebular radius. 

There is no dependence of the magnitude of the ADF upon the excitation energy of 
the UV, optical or IR CEL transition used, indicating that classical nebular temperature 
fluctuations — i.e. in a chemically homogeneous medium — are not the cause of the observed 
abundance discrepancies. Instead, we conclude that the main cause of the discrepancy is en- 
hanced ORL emission from cold ionized gas located in hydrogen-deficient clumps inside 
the main body of the nebulae, as first postulated by Liu et al. (2000) for the high-ADF PN 
NGC 6153. We have developed a new electron temperature diagnostic, based upon the rela- 
tive intensities of the O II 4f-3d A4089 and 3p-3s A4649 recombination transitions. For six 
out of eight PNe for which both transitions are detected, we derive ORL electron temper- 
atures of <300K, very much less than the forbidden-line and H+ Balmer jump temper- 
atures derived for the same nebulae. These results provide direct observational evidence for 
the presence of cold plasma regions within the nebulae, consistent with gas cooled largely by 
infrared fine structure and recombination transitions; at such low temperatures recombination 
transition intensities will be significantly enhanced due to their inverse power-law temperature 
dependence, while UV and optical CELs will be significantly suppressed. 
Key Words: ISM: abundances - planetary nebulae: general 

1 INTRODUCTION 



This is the second of two papers devoted to the study of elemental 



2 KG. Tsamis et al. 



tary nebulae (PNe). In a companion paper, Tsamis et al. (2003a) 
have presented a similar analysis of a number of Galactic and Mag- 
ellanic Cloud H II regions. The main focus of these papers is on the 
problem of the optical recombination-line emission from heavy ele- 
ment ions (e.g. C^^, N^^, O^^) in photoionized nebulae. The main 
manifestation of this problem is the observed discrepancy between 
nebular elemental abundances derived from weak, optical recombi- 
nation hnes (ORLs; such as C ii A4267, N ii A4041, O ii AA4089, 
4650) on the one hand and the much brighter collisionally-excited 
lines (CELs; often collectively referred to as forbidden lines) on 
the other (Kaler 1981; Peimbert, Storey & Torres-Peimbert 1993; 
Liu et al. 1995, 2000, 2001b; Gamett & Dinerstein 2001a; Tsamis 
2002; Tsamis et al. 2003a), with ORLs typically being found to 
yield ionic abundances that are factors of two or more larger than 
those obtained from CELs emitted by the same ions. A closely 
Unked problem involves the observed disparity between the neb- 
ular electron temperatures derived from the traditional [O ill] 
(A4959+A5007)/A4363 CEL ratio and the H I Balmer discontinu- 
ity diagnostic: the latter yields temperatures that are in most cases 
lower than those derived from the [O ill] ratio (Peimbert 1971; Liu 
& Danziger 1993b; Liu et al. 2001b; Tsamis 2002). 

The ORL analysis of the current paper is based upon deep op- 
tical spectra of twelve galactic and three Magellanic Cloud PNe 
that were acquired by Tsamis et al. (2003b; hereafter Paper I). Pa- 
per I describes how the observations were obtained and reduced 
and presents tabulations of observed and dereddened relative in- 
tensities for the detected lines. CoUisionally excited lines (CELs) 
in the spectra were used to derive nebular electron temperatures 
and densities from a variety of diagnostic ratios. They also derived 
CEL-based abtmdances for a range of heavy elements, using stan- 
dard ionization correction factor (icf) techniques to correct for un- 
observed ion stages. 

In the current paper we analyze the ORL data that were pre- 
sented in Paper I. In Section 2 we derive recombination-line ionic 
abundances for a mmiber of carbon, nitrogen and oxygen ions. Sec- 
tion 3 presents a comparison between total C, N and O abundances 
derived from ORLs and from ultraviolet, optical and infrared CELs 
and derives abundance discrepancy factors (ADF's; the ratio of the 
abundances derived for the same ion from ORLs and from CELs) 
for a range of carbon, nitrogen and oxygen ions. In Section 4 we 
investigate how ORL/CEL ADF's correlate with other nebular pa- 
rameters, such as the difference between [O III] forbidden line and 
H I Balmer jump temperatures; the H/3 nebular surface brightness; 
and the nebular absolute radius. Section 5 looks at whether the ob- 
servational evidence provides support for the presence of classi- 
cal Peimbert-type temperature fluctuations within the nebulae, and 
whether the observational evidence points to strong density varia- 
tions within the nebulae. In Section 6 we present evidence for the 
presence of cold plasma (Te < 2000 K) in a number of nebulae in 
our sample, making use of the fact that the strengths of several well 
observed O II and He I recombination lines have sufficiently differ- 
ent temperature dependences for the relative intensities of two O ii 
lines, or two He I lines, to be used as diagnostics of the electron 
temperatures prevailing in their emitting regions. Section 7 sum- 
marizes our conclusions. 

2 RECOMBINATION-LINE ABUNDANCES 

2.1 Carbon ions: C'^+m+, C^+m+ and C*+/H+ 

We have detected recombination lines of carbon from all of the 
PNe, with the exception of LMC N66. C II lines in particular were 



detected from all galactic nebulae, as well as from SMC N87 and 
LMC N141; C III lines were detected from the majority of them, 
excluding NGC3132 and My Cn 18 only. The strongest observed 
C II optical recombination line is the A4267 (V6) 3d-4f transi- 
tion, which was consistently recorded with a high signal-to-noise 
(S/N) ratio in our high resolution deep spectra (see Fig. 1 in Pa- 
per I). We derived C^"'"/H+ abundance ratios from it using the 
recently calculated effective recombination coefficients of Davey, 
Storey & Kisielius (2000), which include both radiative and dielec- 
tronic processes; it should be noted that abundance ratios derived 
using A4267 are insensitive to the adoption of Case A or B. The 
C II A4267 line has been used in the past in abundance analyses of 
galactic PNe (e.g. Kaler 1986, Rola & Stasihska 1994) with con- 
flicting results. It has also been detected from a number of LMC 
and SMC planetary nebulae (Barlow 1987; Meatheringham & Do- 
pita 1991; and from SMC N42 by Vassihadis et al. 1992). No detec- 
tions of this line have been reported however for LMC N141 and 
SMC N87 and carbon ORLs have not been used to date in abun- 
dance studies of Magellanic Cloud PNe. 

Due to the structure of the C^ ion, the configuration of the va- 
lence orbital gives rise to only one atomic term, compared to three 
atomic terms for N"*" and O"*" (e.g. Kuhn 1969; Allen 1973). As a 
result there are fewer C II recombination lines than of O II or N II, 
so they are of greater intensity. As an observational consequence 
of this fact, in several PNe of our sample we have also detected 
C II recombination lines from higher principal quantum numbers, 
originating from states above the 4f ^F° level. In Table 1 we com- 
pare the observed intensities of these lines of high excitation en- 
ergy, normalized such that /(A4267) = 1.00, against the recombi- 
nation theory predictions of Davey et al. (2000). This permits us to 
check whether the 4f '^F°-ng transitions which populate the up- 
per level of 3d^f A4267 can be safely attributed to recombination 
only, or whether unidentified processes contribute as well. This is of 
importance in the light of results from this and previous works (e.g. 
the PN studies of Kaler 1981, 1986; Rola & Stasihska 1994 and 
Liu et al. 1995, 2000 as well as the Hll region analysis of Tsamis 
et al. 2003a), which found that the C^^/H^ abundances derived 
from the A4267 recombination line are often significantly higher 
than those derived from the coUisionally excited C ill] A 1908 line. 
This fact had sometimes been attributed in the past to erroneous re- 
combination coefficients, inaccurate line detections, or blending of 
the A4267 transition with a line from an unknown ionic species. 

In the case of NGC 3242, Table 1 shows that the agreement 
between observations and theory is excellent for all detected lines. 
The 4f ^F°-7g^G A5342 line is blended with a feature identified 
as [Kr IV] A5345.9 (cf Hyung, Aller & Feibelman 1999, their Ta- 
ble 6) and its intensity was retrieved though Gaussian profile fit- 
ting. For NGC 5315 the agreement is very good for both detected 
lines of high excitation. For NGC 5882 the A6462 line is stronger 
by 49 per cent and the A4802 line weaker by 19 per cent than pre- 
dicted, relative to A4267. For IC 4191, as measured on the fixed-slit 
spectrum, the former line is stronger than predicted by 40 per cent, 
while the A5342 line is within 8 per cent of the predicted value. It 
should be noted that apart from A4802, the other high-level C II 
lines are covered in our lower resolution 4.5 A FWHM spectra only 
and the modest discrepancies for NGC 5882 and IC 4191 regarding 
the A6462 line are within the estimated error margins. In the above 
cases the A4802 line was deblended from the N ii A4803.3 line via 
Gaussian line fitting. 

In conclusion, the good consistency among the observed and 
predicted relative intensities of C ii 4f ^F°-ng ^G transitions in this 
PN sample suggests that there is no other mechanism competing 
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with recombination that could contribute to the excitation of the 
3d-4f A4267 line. Thus the high S/N ratio detections of the other 
C II lines from the current PN sample indicate that C ii A4267 is 
a reliable abundance diagnostic. A similar conclusion can be 
drawn from the C II recombination lines found in the spectrum of 
the Orion Nebula (M42). In Table 1 we extend the above compar- 
ison to 4f-ng C II transitions from M 42, whose intensities were 
presented recently by Baldwin et al. (2000). These authors did not 
mark the listed lines as C ii transitions, but left them unidentified 
(cf. their Table 1). However, the measured wavelengths and our 
examination of their relative intensities leave no doubt as to their 
identity. The agreement with theory is excellent in this case as well, 
confirming beyond reasonable doubt the interpretation of recombi- 
nation excitation for these lines. 

The C^"'"/H+ abundance ratios were derived from the A4187 
(V 18) line and from the A4650 (V 1) multiplet, using the effective 
radiative and dielectronic recombination coefficients of Pequignot, 
Petitjean & Boisson (1991) and Nussbaumer & Storey (1984), re- 
spectively; they are insensitive to the assiraiption of Case A or B. 
The C III (V 16) triplet at A4069 was also detected from NGC 2022, 
NGC3242, NGC 5882 and NGC 68 18. This multiplet is primarily 
excited by radiative recombination and is usually seriously blended 
with [S II] and O ii V 10 lines in medium resolution spectra; its 
overall intensity was retrieved via multiple Gaussian fitting. As was 
noted by Liu (1998), the observed relative intensities of C ill V 16, 
V 18 and V 1 multiplets are not in accord with theoretical predic- 
tions; in an analysis of the spectrum of NGC 4361 he found that the 
(2;3+/jj+ abundance ratio derived from these three C iii multiplets 
spans a range of 0.4 dex, which he attributed to probable uncertain- 
ties in the effective recombination coefficients. The four PNe for 
which all V 16, V 18 and V 1 multiplets have been detected offer 
a further testing ground for the reliability of current ORL C III re- 
combination coefficients. An instructive case is that of the high ex- 
citation nebula NGC 2022. As is apparent from our high resolution 
1.5 A FWHM spectra of this object, the auroral [S II] AA4068, 4075 
lines are virtually absent and the observed intensity of C III V 16 
is most accurately recorded, since it is less affected by blending ef- 
fects. In this nebula, the observed C iii AA4069, 4187 and A4650 
multiplets have intensity ratios of 1.80 : 0.34 : 1.00, versus the the- 
oretical ratios of 0.59 : 0.21 : 1.00 calculated using the effective 
radiative and dielectronic recombination coefficients of Pequignot 
et al. (1991) and Nussbairaier & Storey (1984), respectively. In the 
case of NGC 5882, the observed multiplet intensities show ratios 
of 2.38 : 0.22 : 1.00. Similarly for IC4191, where C III A4069 is 
not detected, the observed intensities of A4187 (V 18) and A4650 
(V 1) have a ratio of 0.20 : 1.00. We see that the theoretical predic- 
tions for the relative intensities of A4187 and A4650 appear to be 
more secure than that for A4069. It would seem that the effective 
recombination coefficient for the A4069 V 16 multiplet has been 
underestimated by a factor of ~2-4, or else that it is blended with 
an unknown line from an ion of similar excitation. 

As a result of the above discussion, in our ORL abundance 
analysis for C^+ we have discarded values derived from the A4069 
triplet; instead, we have adopted C^"'"/H+ abundance ratios as de- 
rived from an intensity- weighted mean of the C ill A4187 and 
A4650 lines. 

Regarding C*"*", we have derived C*"'"/H"'" fractions for a num- 
ber of galactic PNe from the C IV A4658 line, using the Case 
A effective recombination coefficients of Pequignot et al. (1991). 
The ionic and total carbon abundances derived from recombination 
lines are presented in Tables 2 and 3 for the Galactic and Magel- 



Table 3. Carbon and oxygen recombination-line abundances for Magellanic 
Cloud PNe. 





Cl^Xf^ TVTOT 

oMC JNov 




IjMC JM41 


/(A4267) 


0.672 


* 


0.683 


10'' X C2+/H+ 


6.73 


* 


6.79 


/(A4187) 


0.108 


* 


0.0522 


W X C3+/H+ 


1.46 


* 


0.733 


7(A4650) 


0.141 


* 


0.181 


10'' X C^+/H+ 


0.413 


* 


0.538 


Adopted 








10^ X C3+/H+ 


6.75 


* 


5.82 


/(A4658) 


* 


* 


* 


10* X C*+/H+ 


* 


* 


* 


icf(C) 


1.020 


* 


1.023 


10^ X C/H 


7.55 




7.54 


/(A4069)|O2+ 




* 


0.188|7.27 


/{A4072) 02+ 




* 


0.124 5.15 


/(A4075) o2+ 






0.254 7.30 


/(A4638) 0^+ 


.0832|7.73 


* 


.0255 2.38 


7(A4641) 


.0831 3.06 


* 


0.120 4.45 


/(A4649) 0^+ 


.0781 1.51 


0.227|3.98 


0.127 2.47 


/(A4650) 02+ 


.0781 7.26 


0.227|19.1 


.08367.81 


/(A4661) 0^+ 


.0243 1.76 


* 


.0518 3.79 


/(A4676) 0^+ 


.0248 2.15 


* 


.0465 4.04 


Adopted 








lO"* X 02+/H+ 


2.95 


8.53 


4.96 


CEL 0+/02+ 


0.018 


0.076 


0.023 


icf(0) 


1.000 


3.770 


1.000 


10" X O/H 


3.00 


34.6 


5.08 



lanic Cloud PNe, respectively (see Appendix A for a discussion of 
the adopted icf scheme in the derivation of total C abundances). 

2.2 Nitrogen ions: N^+/H+, N^+/H+ and N^+/H+ 

We have detected recombination lines from all galactic PNe and 
from up to three ionization stages of nitrogen. Nitrogen ORLs 
were not detected from any of the three Cloud nebulae. Both 
3s-3p as well as 3d-4f transitions were recorded. The derived 
recombination-line N^+ , N'^+ , N''+ and total N abundances are pre- 
sented in Table 4 (see also Appendix A). 

The N^^/H^ fractions were derived using effective recom- 
bination coefficients from Escalante & Victor (1990), assimring 
Case A for singlets and Case B for triplets. The adopted N^+/H+ 
abundances are derived for each nebula by averaging the val- 
ues obtained from each N II line, weighted according to the pre- 
dicted relative intensity of the transition. The A4379 (V 18) N III 
recombination line was detected from all galactic PNe, except 
NGC 5315 and My Cn 18. We derived N^+/H+ abundances from 
it using the effective radiative and dielectronic recombination co- 
efficients of Pequignot et al. (1991) and Nussbaumer & Storey 
(1984) respectively. The N IV 5g-6h A4606 line has also been 
detected from the high excitation nebulae NGC 2022, NGC 2440, 
NGC 3918, NGC 6302 and NGC 6818; it was used for the deriva- 
tion of N''^/H^ abundances employing effective recombination co- 
efficients as in the case of N ill A4379. 

2.3 Oxygen ions: 0^+m+, 0^+m+ and 0''+/H+ 

Our observations provide us with one the most extensive records 
of O II recombination spectra in gaseous nebulae thus far, obtained 



4 Y. G. Tsamis et al. 



Table 1. High-excitation C II recombination line relative intensities. 



Ao(A) 


Trans. 


NGC 3242 


NGC 53 15 


NGC 5882 


IC4191 


Orion (M42) 


Theory 






-^obs 


-^obs 


-^obs 


-^obs 


-^obs 


^pred 


4267.15 


3d-4f 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


6151.43 


4d-6f 


0.030 


* 




*|0.055" 


0.038 


0.040 


6461.95 


4t-6g 


0.102 


0.095 


0.153 


0.166|0.144 


0.087 


0.103 


5342.38 


4f-7g 


0.065 


* 




*|0.056 


0.049 


0.053 


4802.23 


4f-8g 


0.031 


0.020 


0.025 




0.034 


0.031 



° The value before the dash is for the entire nebula; those after are from a fixed-slit observation. 



Table 2. Recombination line carbon abundances for Galactic PNe. 





N(iC 3242 




N(iC 5315 


N(iC 3918 


NGC 2022 




/(A4267) 


0.620 


0.399 


0.706 


0.497 


0.820 


0.449 


1 a4 . - /^24- /it4- 
10 X C^/H^ 


6.15 


3.77 


6.60 


5.02 


8.83 


4.69 


/(A4069) 


0.755 


0.412 


* 


* 


1.75 


0.957 


10-* X C3+/H+ 


3.64 


2.11 


* 


* 


7.49 


3.97 


7(A4187) 


0.204 


0.0376 


0.0107 


0.152 


0.332 


0.118 


10 X C^+/H+ 


2.83 


0.552 


0.161 


2.03 


4.07 


1.41 


7(A4650) 


0.590 


0.173 


* 


0.418 


0.974 


0.389 


10* X C3+/H+ 


1.75 


0.522 


* 


1.22 


2.77 


1.10 


Adopted 














10* X C3+/H+ 


2.03 


0.527 


0.161 


1.44 


3.10 


1.37 


/(A4658) 


0.0776 




* 


0.109 


1.03 


0.180 


10* X C*+/H+ 


0.188 




* 


0.262 


2.62 


0.447 


icf(C) 


1.010 


1.030 


1.078 


1.116 


1.020 


1.185 


10* X C/H 


8.45 


4.43 


7.29 


7.50 


14.84 


7.71 




NGC 3132 


NGC 2440 


NGC 6302 


IC4406 


IC4191 


MyCn 18 


7(A4267) 


0.697 


0.403 


0.163 


0.805 


0.546 


0.252 


10* X C2+/H+ 


6.60 


4.47 


2.00 


7.72 


5.16 


3.80 


7(A4187) 


* 


0.222 


* 


0.0912 


0.0711 


* 


10* X C3+/H+ 


* 


2.58 


* 


1.32 


1.04 


* 


7(A4650) 


* 


0.337 


0.0754" 


0.238 


0.359 


* 


10* X C3+/H+ 


* 


1.06 


* 


0.715 


1.08 


* 


Adopted 














10* X C3+/H+ 


* 


1.66 


* 


0.883 


1.07 


* 


/(A4658) 




0.722 




* 


* 




10* X C*+/H+ 




1.82 




* 


* 




icf(C) 


1.930 


1.284 


2.725 


1.310 


1.040 


1.589 


10* X C/H 


12.74 


10.21 


5.45 


11.27 


6.48 


6.04 



« From C III A4647. 



for PNe possessing a wide range of physical conditions. For the first 
time also to our knowledge, O 11 lines have been detected and mea- 
sured from Magellanic Cloud planetary nebulae (SMC N87, LMC 
N66, and LMC N141). This allowed us to obtain accurate ORL 
Q2-i-yjj+ abundances for all our nebulae and permitted a compre- 
hensive investigation on the occurrence of the discrepancy between 
recombination-line abundances and those derived in the usual 
manner from forbidden lines of O^"*". Furthermore, having such an 
extensive inventory of O ii line intensities, we are able to perform 
a thorough comparison with current theoretical predictions of the 
O II recombination spectrum (see Appendix B). 

In Table 5 we present ORL O'^^/H^ ionic ratios for the com- 
plete sample (for the Magellanic Cloud objects, these ratios and to- 
tal O abtmdances are listed in Table 3). Effective recombination co- 



efficients are taken from Storey (1994) for 3s-3p transitions (under 
LS-coupling) and from Liu et al. (1995, hereafter LSBC) for 3p- 
3d and 3d-4f transitions (under intermediate coupling), assuming 
Case A for doublet and Case B for quartet lines. For several nebu- 
lae, we have also detected lines arising from doubly excited spec- 
tral terms, such as multiplet V 15 3s' ^D-3p' ^F° at 4590 A (cf. Pa- 
per I). The excitation of this multiple! is dominated by dielectronic 
recombination, but we have not derived abundance ratios from it 
since the existing recombination coefficients are not of the desired 
accuracy.^ For each PN, the mean 0^"'"/H+ fractions derived by av- 

^ Gamett & Dinerstein (2001a) noted that for a sample of about 12 PNe, 
abundances derived firom multiplet V 15 of O ii were anomalously 
high compared to those derived from other O II multiplets, owing either to 
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Table 4. Recombination-line N-^+/H+ abundances for Galactic PNe. 









N2+ 




N2+ 




N2+ 




N2 + 


Ao 


Mult. 


-^obs 


H+ 


-^obs 


H+ 


-^obs 


H+ 


-^obs 


H+ 


(A) 






(10-4) 




(10-4) 




(10-4) 




(10-4) 






NGC3242 


NGC5882 


NGC5315 


NGC3918 


V3 3s3p°- 


-3p3D 


















5666.63 


V3 


.0321 


2.34 






* 


* 


* 


* 


5676.02 


V3 


.0134 


2.20 


.0189 


3.10 


* 


* 


* 


* 


5679.56 


V3 


.0428 


1.68 


.0604 


2.36 


* 


* 


* 


* 


V5 3s3p°. 


-3p3p 


















4601.48 


V5 


.0050 


1.19 


* 


* 


* 


* 


* 


* 


4607.16 


V5 


* 


* 


* 


* 


* 


* 


* 


* 


4621.39 


V5 


* 


* 


* 


* 


* 


* 


* 


* 


4630.54 


V5 


.0179 


1.43 


* 


* 


* 


* 


* 


* 


3d-4f 




















4041.31 


V39b 


.0128 


0.88 


.0273 


1.79 


.0573 


3.71 


.0134 


0.94 


4237.05 


V48b 


.0173 


1.71 


.0139 


1.38 


.0166 


1.64 


.0163 


1.65 


4241.78 


V48a 










.0422 


4.16 


.0163 


1.75 


4530.41 


V58b 






.0180 


1.75 


=1: 


si; 


.0148 


1.55 


4678.14 


V61b 










.0312 


4.80 




Hi 


Sum 




0.141 


1.71 


0.138 


2.24 


0.147 


3.43 


.0607 


1.42 






NGC6302 


IC4191 


NGC6818 


NGC2022 


V3 3s3p°. 


-3p3D 


















5666.63 


V3 


* 


* 


.0921 


6.90 






* 


* 


5676.02 


V3 


* 


* 


.0400 


6.58 


.0322 


5.29 


.0288 


4.73 


5679.56 


V3 


* 


* 


.1282 


5.02 


.1029 


4.03 


.0921 


3.61 


V5 3s ^P°- 


-3p^P 


















4630.54 


V5 


.0743 


5.55 


.0549 


4.50 


* 


* 


* 


* 


V 28 3s P 


-3p D 


















5927.81 


V28 


.0205 


8.26 






* 


* 


* 


* 


5931.78 


V28 


.0361 


6.47 


.0460 


5.73 




* 


* 


* 


3d^f 




















4041.31 


V39b 






.0639 


4.26 






.0580 


4.27 


Sum 




0.131 


5.64 


0.425 


5.50 


0.135 


4.27 


0.179 


4.38 






NGC3132 


NGC2440 


IC4406 


MyCnlS 


V3 3s3p°- 


-3p '^D 


















5676.02 


V3 


.0249 


4.10 














5679.56 


V3 


.0799 


3.13 


* 




* 


* 


* 


* 


V5 3s3p°. 


-3p3p 


















4607.16 


V5 


* 


* 


* 


* 


* 


* 


.0842 


26.7 


4621.39 


V5 


* 


* 


* 


* 


* 


* 


.0731 


23.3 


4630.54 


V5 


.0411 


3.40 


* 


* 


.2318 


1.90 


.2099 


17.9 


3d-4f 




















4041.31 


V39b 






.0618 


4.28 










4237.05 


V48b 






.0495 


4.90 










4241.78 


V48a 






.0558 


5.89 






* 


* 


4530.41 


V58b 


* 


* 


.0415 


4.28 


* 


* 


* 


* 


Sum 




0.146 


3.40 


0.209 


4.77 


0.232 


1.90 


0.283 


20.4 



eraging the values from all 3-3 multiplets and the co-added 3d-4f 
transitions have been adopted as the recombination line values in 
the subsequent discussion. The O^^/H^ values listed for the 3d-4f 
lines in Table 6 were obtained by summing all the intensities and 
dividing by the sum of all the recombination coefficients, allowing 



an underestimated dielectronic recombination coefficient for this multiplet, 
or to additional contributions, e.g. high-temperature dielectronic recombi- 
nation originating from very hot regions, which might also contribute to the 
strengths of other O II Unes. Liu et al. (2001b) looked at this issue in de- 
tail and concluded that the discrepancy was most likely due to the lower 
accuracy of the recombination coefficients available for the V 15 multiplet. 



for weak unobserved or blended multiplet components, as listed in 
Table 4(a)ofLSBC. 

In Table 7 we present a simmiary of the oxygen 
recombination-line ionic and total abundances, including the 
q3-i-^j_j-i- fi-aQ(ions for several PNe derived by Liu & Danziger 
(1993a) from the O 111 V8 multiplet at 3265 A; we also hst 
Q'^+/ll+ fractions derived from the O IV A4632 line which was 
detected from the high excitation nebulae NGC3918, NGC2022, 
NGC 6818, NGC 2440 and NGC 6302. For the latter calculation we 
used the Case A dielectronic recombination coefficients of Nuss- 
baumer & Storey (1984) and the radiative recombination coeffi- 
cients of Pequignot et al. (1991). For a discussion on the derivation 
of total O abundances the reader is referred to Appendix A. 
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Table 5. Overall recombination-line nitrogen abundances for Galactic PNe. 





NGC3242 


NGC 5882 


NGC 5315 


NGC 3918 


NGC 2022 


NGC 6818 


10* X N2+/H+ 


1.71 


2.24 


3.43 


1.42 


4.38 


4.27 


/(A4379) 


0.158 


0.136 


* 


0.146 


0.472 


0.203 


10" X N3+/H+ 


0.685 


0.575 


* 


0.641 


2.080 


0.887 


/(A4606) 


* 


* 


* 


.0660 


0.112 


.0460 


10* X N4+/H+ 


* 


* 


* 


0.207 


0.368 


0.146 


CEL N+/N2+ 


0.008 


0.028 


0.163 


0.203 


0.013 


0.241 


10"* X N/H 


2.42 


2.88 


5.93 


2.56 


6.88 


6.33 




NGC3132 


NGC 2440 


NGC 6302 


IC4406 


IC4191 


MyCnlS 


10"* X N2+/H+ 


3.40 


4.77 


5.64 


1.90 


5.00 


20.35 


/(A4379) 


.0571 


0.324 


0.182 


0.172 


0.270 




10* X N3+/H+ 


0.247 


1.430 


0.807 


0.731 


1.150 


* 


/(A4606) 


* 


0.475 


0.668 


* 


* 


* 


10* X N*+/H+ 


* 


1.70 


2.28 




* 


* 


CEL N+/N^+ 


1.296 


0.475 


0.435 


0.536 


0.393 


* 


10* X N/H 


8.05 


10.17 


11.18 


3.65 


8.12 


>20.35 



Table 7. Recombination-line oxygen abundances for Galactic PNe 





NGC 3242 


NGC 5882 


NGC 5315 


NGC 3918 


NGC 2022 


NGC 6818 


10* X 02+/H+ 


6.28 


9.70 


8.57 


5.36 


13.01 


7.33 


10* X 03+/H+ 


2.05" 


* 


* 




3.62" 


2.37" 


/(A4632) 


* 


* 


* 


0.150 


0.334 


0.200 


10* X 0*+/H+ 


* 


* 


* 


0.474 


1.100 


0.641 


CEL 0+/02+ 


0.0091 


0.029 


0.098 


0.086 


0.020 


0.090 


icf(0) 


1.000 


1.010 


1.332 


1.298 


1.000 


1.000 


10* X O/H 


8.39 


10.08 


12.53 


12.32 


17.99 


11.00 




NGC 3132 


NGC 2440 


NGC 6302 


IC4406 


IC4191 


MyCnl8 


10* X 02+/H+ 


8.15 


5.23 


3.28 


7.06 


12.92 


6.43 


10* X 03+/H+ 


* 


2.56" 


* 


* 


* 


* 


/(A4632) 


* 


0.300 


0.228 


* 


* 


* 


10* X 0*+/H+ 


* 


1.01 


0.784 


* 


* 


* 


CEL 0+/02+ 


0.894 


0.285 


0.092 


0.471 


0.041 


0.590 


icf(0) 


1.020 


1.000 


1.672 


1.040 


1.070 


1.003 


10* X O/H 


15.74 


10.29 


7.30 


10.80 


14.39 


10.24 



" From the Liu & Danziger (1993a) analysis of O III V8 A3265. 



3 COMPARISON OF ORL AND CEL ABUNDANCES 
3.1 Total C, N and O abundances 

The total elemental abundances of C, N and O obtained for the 
current PN sample from CELs (Paper I) and from ORLs are pre- 
sented in Table 8, on the usual log iV(H) = 12.0 scale, where they 
are compared with average abundances for Galactic PNe (for both 
Type I and non-Type I nebulae) derived by Kingsburgh & Barlow 
(1994) and the solar photospheric abundances of Grevesse, Noels 
& Sauval (1996) and Allende Prieto et al. (2001, 2002). 



3.2 Abundance grids and ionic ADF's 

Comparisons between the ionic abundances obtained from the UV, 
optical and IR colUsionally excited Unes (Paper I) and those from 
ORLs (Tables 2-7) are presented in Fig. 1 bringing together lUE, 



optical, and ISO data for 9 galactic PNe, and WE and optical data 
for the 3 Cloud PNe (LWS spectra are not available for NGC 2440 
and My CN 1 8, while lUE spectra are not available for IC 4191 and 
MyCN 18). Note that the plotted 0^+/H+ fractions derived from 
optical CELs are those from the [O III] AA4959, 5007 lines only. 

In Paper I we used the ISO LWS far-IR fine structure (FS) lines 
presented by Liu et al. (2001a) to derive N2+/H+ and 0^+/H+ 
abundances for 8 PNe in our sample. In the meantime we re- 
trieved archived LWS spectra of NGC 2022 (TDT 69201604) and 
NGC 6818 (TDT 34301005) and measured the fluxes of the [N III] 
57-/im and [O ill] 52-, 88-/xm FS lines. We scaled these to /(H/3) = 
100 using the total nebular H/3 fluxes dereddened by the c(H/3)''"*° 
constants of Table 5, Paper I. The respective /(52-/im), 7(57-/im), 
and J(88-/xm) intensities are 9L7, 21.5, and 42.2 for NGC 2022; 
and 164.6, 45.5 and 69.7 for NGC 6818. From the [O iii] 52- 
/^m/88-/^m line ratio we derived electron densities of 800 cm~^ for 



A deep suryey of heavy element lines in PNe - II 

Table 6. Recombination-line 0'^+/H+ abundances for galactic Table 6. -continued 



PNe. 















Ao 


Mult. 




02+ 




02+ 


Ao 


Mult. 


-^obs 








-^obs 


H+ 


-^obs 


H+ 


H+ 


-^obs 


H+ 


(A) 






(10-4) 




(10-4; 








(10-4) 




(10-4) 




























NGC3242 


NGC5882 






NGC3242 


NGC5882 


4489.49 


V86b 


.0055 


7.86 


* 


* 


4638.86 


VI 


.1692 


15.9: 


.1690 


16.2: 


4491.23 


V86a 


.0235 


16.0 


.0203 


13.3 


4641.81 


VI 


.2152 


8.00 


.2522 


9.59 


4609.44 


V92a 


.0439 


7.26 


.0609 


9.64 


4649.13 


VI 


.2231 


4.36 


.3447 


6.89 


4669.27 


V89b 


.0119 


27.6: 




* 


4650.84 


VI 


.0979 


9.18 


.1141 


10.9 


3d-^f 




0.493 


8.40 


0.569 


9.67 


4661.63 


VI 


.0870 


6.38 


.1188 


8.92 


Adopted 






6.28 




9.70 


4673.73 


VI 


.0304 


13.1: 


.0188 


9.10 






NGC5315 


NGC3918 


4676.24 


VI 


.0574 


5.01 


.0906 


8.10 


4638.86 


VI 


.1740 


16.5: 


.2200 


20.4: 


V13s4p-3p4D° 




0.681 


6.02 


0.939 


8.29 


4641.81 


VI 


.2367 


9.86 


.2168 


7.96 


4317.14 


V2 


.0454 


5.39 


.0638 


8.43 


4649.13 


VI 


.3574 


7.08 


.2108 


4.07 


4319.63 


V2 


.0227 


2.74 


.0333 


4.08 


4650.84 


VI 


.0900 


8.56 


.0415 


3.84 


4325.76 


V2 


.03 1 1 


20.3: 






4661.63 


VI 


.0897 


6.67 


.0743 


5.39 


4345.56 


V2 


.0881 


11.1: 


.1100 


14.6: 


4673.73 


VI 


.0137 


6.56 


.0310 


14.5: 


4349.43 


V2 


.0712 


3.71 


.0911 


4.82 


4676.24 


VI 


.0789 


6.99 


.0656 


5.66 


4366.89 


V2 


.0909 


10.6: 


* 


* 


V13s4p-3p4D° 




0.866 


7.58 


0.609 


5.28 


V2 3s*P-3p4p° 




0.139 


3.88 


0.188 


5.43 


4317.14 


V2 


.0694 


8.23 


.0417 


4.93 


4414.90 


V5 


.0349 


6.20 


.0439 


8.51 


4319.63 


V2 


.0491 


5.94 


.0170 


2.04 


4416.97 


V5 


.0331 


10.6 


.0399 


13.9 


4325.76 


V2 


.0191 


11.8: 


.0141 


9.2: 


4452.37 


V5 


.0179 


28.8: 


* 


* 


4345.56 


V2 


.0635 


8.34 






V5 3s2p-3p2D° 




0.068 


7.77 


0.084 


10.4 


4349.43 


V2 


.0744 


9.77 


.0638 


3.31 


4069.89 


VIO 


.1828 


7.08 


.1361 


5.27 


4366.89 


V2 


.0589 


6.87 


* 


* 


4072.16 


VIO 


.2153 


8.95 


.2385 


9.91 


V2 3s4p-3p4p° 




0.315 


6.06 


0.123 


3.39 


4075.86 


VIO 


.1664 


4.79 


.2754 


7.92 


4414.90 


V5 


.0721 


13.0 


.0518 


8.88 


4078.84 


VIO 


.0213 


5.81 


* 


* 


4416.97 


V5 


.0620 


20.2 


.0290 


8.96 


4085.11 


VIO 


.0351 


7.81 


.0290 


6.45 


4452.37 


V5 


.0146 


29.8: 


.0653 


100.: 


4092.93 


VIO 


.0240 


7.32 


.0212 


6.45 


V5 3s2p-3p2D° 




0.134 


15.55 


0.081 


8.91 


V10 3p''D°-3d*F 




0.645 


6.71 


0.700 


7.58 


4069.62 


VIO 


.4263 


16.5: 






4121.46 


V19 


.0283 


9.81 


.0289 


10.31 


4072.16 


VIO 


.1879 


7.78 






4129.32 


V19 


.0517 


79.2: 






4075.86 


VIO 


.1412 


4.05 




* 


4132.80 


V19 


.0210 


3.87 


.0535 


9.77 


4085.11 


VIO 


.0437 


9.69 


* 


* 


4153.30 


V19 


.0304 


3.91 


.0955 


12.2 


V10 3p4D°-3d4F 




0.329 


5.58 


* 


* 


4156.53 


V19 


.0373 


30.1: 


.0728 


58.3: 


4121.46 


V19 


* 


* 


.0920 


33.1: 


4169.22 


V19 


.0126 


4.76 


.0279 


10.5 


4132.80 


V19 


* 


* 


.0180 


3.32 


V19 3p4p°-3d4p 




0.091 


4.90 


0.206 


10.9 




viy 


* 


* 


.0260 


3.35 


4110.78 


V20 


* 


* 


.0359 


14.9 


4156.53 


V19 


* 


* 


.0370 


29.9: 


4119.22 


V20 


* 


* 


.0919 


10.3 


4169.22 


V19 


* 


* 


.0045 


1.72 


V20 3p4p°-3d4D 




* 


* 


0.128 


11.3 


V 19 3p F -3d P 




* 


* 


0.049 


3.07 


4890.86 


V28 


* 


* 


.0086 


7.39 


4083.90 


V48b 


.0241 


7.80 






4906.83 


V28 






.0362 


14.5 


4087.15 


V48c 


.0319 


10.9 






4924.53 


V28 






.0656 


15.4 


4Uoy.2y 


V48a 


.0899 


8.30 


* 


* 


V28 3p4p°-3d4D 




* 


* 


0.110 


13,9 


4275.55 


V67a 


.1005 


8.11 


* 


* 


4083.90 


V48b 


.0211 


6.88 


.0210 


6.58 


4303.82 


V53a 


.0331 


6.50 


* 


* 


4087.15 


V48c 


.0204 


7.04 


.0359 


11.8 


4609.44 


V92a 


.0476 


7.79 


* 


* 


4089.29 


V48a 


1113 


10 4 


1452 




3d-^f 




0.332 


8.08 


* 


* 


4275.55 


V67a 


.0310 


4.47 


.0579 


7.99 


Adopted 






8.57 




5.36 


4276.75 


V67b 


.0276 


8.36 


.0236 


7.65 






















4277.43 


V67c 


.0172 


6.90 


.0265 


10.2 


4638 86 


VI 


.1051 


9.41: 


.2160 


19.8: 


4282.96 


V67c 


.0086 


5.24 




* 






4641.81 


VI 


.2538 


9.01: 


.4956 


18.0: 


4283.73 


V67c 


.0095 


9.26 










4649.13 


VI 


.3674 


6.86 


.1687 


3.23 


4285.69 


V78b 


.0120 


5.90 










4650.84 


VI 


.0892 


8.00 


.0262 


2.50 


4288.82 


V53c 


.0055 


5.11 


* 


* 


.0221 


12.2 


4661.63 


VI 


.1093 


7.67 


.0851 


6.11 


4291.25 


V55 


.0162 


9.27 


4676.24 


VI 


.0772 


6.44 




* 


4292.21 


V78c 


.0086 


8.67 


.0078 


7.56 


V13s4p-3p4D° 


1.00 


7.69 


0.280 


3.65 


4294.78 


V53b 


.0314 


10.2 


.0156 


4.88 




V5 


4414.90 


.1247 


19.3 


.0594 


9.90 


4303.83 


V53a 


.0582 


11.6 


.0573 


10.9 


V5 3s2p-3p2D° 


0.125 


19.3 


0.059 


9.90 


4307.23 


V53b 


.0087 


7.71 


* 


* 


4069.78 


VIO 


* 


* 


0.132 


5.11 


4313.44 


V78a 


.0086 


6.59 


* 


* 


.3164 


13.1 


.0150 


13.2 


4072.16 


VIO 


0.175 


7.26 


4353.59 


V76c 


.0098 


8.98 


29.3: 


4075.86 


VIO 


.2949 


8.43 


0.438 


12.6 


4357.25 


V63a 


.0182 










4085.11 


VIO 




* 


.0722: 


15.8: 


4366.53 


V75a 


* 


* 


.0605 


7.13 








.0260 


23.7: 


V10 3p4D°-3d4F 




0.611 


10.3 


0.745 


6.41 


4466.42 


V86b 


.0305 


26.7: 




4477.90 


V88 


.0145 


15.7 


* 


* 
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Table 6. -continued 



Ao 




Mult. 




02+ 




02+ 




-^obs 


H+ 


-^obs 


H+ 










(10-4) 




(10-4 








NGC2022 


NGC 6818 


4089.29 




V48a 


.2020 


16.9 


.0984 


9.37 


4276.75 




V67b 


.1856 


13.8 






4609.44 




V92a 


.0578 


10.0 






3d^f 






0.416 


14.8 


0.098 


9.37 


Adopted 








13.0 




7.33 








NGC3132 


IC4406 


4638.86 




VI 


.1043 


10.0 


.1817 


17.3: 


4641.81 




VI 


.1706 


6.49 


.2627 


9.94 


4649.13 




VI 


.1649 


3.30 


.2585 


5.14 


4650.84 




VI 


.1226 


11.8 


.0959 


9.16 


4661.63 




VI 


.1157 


8.70 


.0988 


7.39 


4676.24 




VI 


.0762 


6.82 


.0481 


4.28 


V13s*P-3p 


4j)0 




0.754 


7.06 


0.764 


6.84 


4069.62 




VIO 


.4086 


15.8: 


.4106 


15.9: 


4072.16 




VIO 


.1484 


6.17 


.3776 


15.7: 


4075.86 




VIO 


.3592 


10.3 


.2530 


7.28 


4085.11 




VIO 


.0433 


9.64 


* 


* 


V 10 3p ''D°- 


-3d*F 




0.551 


8.70 


0.253 


7.28 


Adopted 








8.15 




7.06 








MyCnlS 


NGC 6302 


4649.13 




VI 


.2294 


4.66 


.154 


2.70 


4650.84 




VI 


.1531 


14.9 


.0609 


5.13 


4661.63 




VI 


* 




.0353 


2.33 


Adopted 








6.43 




3.28 










IC4191 










entire nebula 


fixed slit 


4638.86 




VI 


.1912 


18.0: 


.2426 


22.8: 


4641.81 




VI 


.3218 


12.0 


.2020 


12.2 


4649.13 




VI 


.5439 


10.7 


.5757 


11.9 


4650.84 




VI 


.0997 


9.38 


.1168 


14.5 


4661.63 




VI 


.1539 


11.3 


.1442 


10.6 


4676.24 




VI 


.1226 


10.8 


.1084 


9.29 


V13s*P-3p 


4j)o 




1.24 


10.9 


1.15 


10.1 


4317.14 




V2 


* 


* 


.0676 


8.77 


4319.63 




V2 


* 


* 


.0551 


6.62 


4345.56 




V2 


* 


* 


.1210 


15.2 


4349.43 




V2 


* 


* 


.2324 


12.1 


4366.89 




V2 


* 


* 


.1064 


12.4 


V2 3s^P-3p 


4po 




* 


* 


0.582 


11.1 


4414.90 




V5 


.0757 


13.2 


.0546 


9.49 


4416.97 




V5 


.0531 


16.6 


.0444 


13.9 


V5 3s2p-3p 


2D° 




0.129 


14.4 


0.099 


11.1 


4069.62 




VIO 


.7056 


27.2: 


.5998 


23.2: 


4072.16 




VIO 


.3654 


15.2 


.3565 


14.8 


4075.86 




VIO 


.2883 


8.28 


.3452 


9.91 


4085.11 




VIO 


.0503 


11.0 


.0266 


5.81 


V10 3p''D°- 


Jd^F 




0.704 


11.1 


0.728 


11.5 


4083.90 




V48b 


.0418 


13.8 


.0701 


23.2 


4087.15 




V48c 


.0580 


18.6 


.0685 


21.9 


4089.29 




V48a 


.1475 


13.9 


.2059 


19.1 


4275.55 




V67a 


.2096 


16.2 


.2434 


18.8 


4282.96 




V67c 


.0587 


12.6 


* 


* 


4303.82 




V53a 


.0930 


18.1 


.0684 


13.3 


4466.42 




V86b 


* 


* 


.0387: 


35.9: 


4609.44 




V92a 


.0852 


14.3 


.0603 


10.1 


3d-^f 






0.694 


15.2 


0.717 


17.6 


Adopted 








12.9 




12.3 



Table 6. -continued 











02+ 






Ao 




Mult. 


-^obs 


H+ 






(A) 






(10-4) 












NGC 2440 






4638.86 




VI 


.2840: 


24.0: 






4641.81 




VI 


.2051: 


7.14: 






4649.13 




VI 


.1266 


2.32 






4650.84 




VI 


.0510 


4.48 






4661.63 




VI 


.0346 


2.38 






4676.24 




VI 


.0710: 


5.81: 






V13s4p-3p 


4jjo 




0.212 


2.63 






4414.90 




V5 


.0472 


6.91 






4416.97 




V5 


.0380 


10.0 






V5 3s2p-3p 






0.085 


8.02 






4153.30 




V19 


.0260 


3.35 






VWSp"?"- 






0.026 


3.35 






Adopted 








5.23 






Table 8. Elemental C, N, O abundances by number, derived from CELs 


and ORLs, in units where log 


Af(H) = 


12.0." 








PN 




C 




N 


O 






ORL 


CEL 


ORL 


CEL 


ORL 


CEL 


NGC 2022 


9.17 


8.33 


8.84 


7.46 


9.26 


8.66 


NGC 2440© 


9.01 


8.37 


9.01 


8.16 


9.01 


8.39 


NGC 3132 


9.11 


8.50 


8.91 


8.37 


9.20 


8.82 


NGC 3242 


8.93 


8.14 


8.38 


7.53 


8.92 


8.52 


NGC 3918 


8.88 


8.64 


8.41 


8.02 


9.09 


8.86 


NGC 5315 


8.86 


8.33 


8.77 


8.52 


9.10 


8.79 


NGC 5882 


8.65 


8.18 


8.46 


8.18 


9.00 


8.69 


NGC 6302(1) 


8.74 


7.89 


9.05 


8.52 


8.86 


8.40 


NGC 6818 


8.89 


8.41 


8.80 


7.74 


9.04 


8.71 


IC4191 


8.81 




8.91 


7.59 


9.16 


8.78 


IC4406 


9.05 


8.56 


8.56 


8.32 


9.03 


8.76 


MyCnlS 


8.78 


* 


>9.31 


8.34 


9.01 


8.75 


SMCN87 


8.88 


8.58 




7.55 


8.48 


8.03 


LMCN66 


* 


7.52 




7.99 


9.54 


8.50 


LMCN141 


8.88 


8.30 




7.95 


8.71 


8.29 


KB94 non-Type l'' 




8.81 




8.14 




8.69 


KB94 Type I*" 


* 


8.48 


* 


8.72 


* 


8.65 


Solar= 


* 


8.39 


* 


7.97 


* 


8.69 



" All values are for the entire nebulae, except NGC 6302, NGC 6818, 
MyCnlS and IC4406 where values are from fixed-slit observations; 

CEL results are from Paper I. 

From Kingsburgh & Barlow (1994). 

Solar photospheric abundances from Grevesse et al. (1996), except 
O and C which are from Allende Prieto et al. (2001, 2002) respectively. 

NGC 2022 and 950 cm"'^ for NGC 6818. From the measured FS 
line intensities and adopting the A'e's from the [O III] FS line ra- 
tio we derived: N2+/H+ = 2.75x10-^ and 0^+/H+ = l.OVx IQ-* 
for NGC 2022; and 6.91x10"-^ and 2.20x10"* respectively for 
NGC 6818. On the other hand adopting the mean of the A^c's from 
the [CI III] and [Ar IV] optical line ratios (Paper 1) we derived: 
N^+/H+ = 3.94x 10-^ and 0^+/H+ = 1.45 x IQ-* for NGC 2022; 
and 1.15x10-4 and 3.41x10-4 respectively for NGC 6818. The 
latter values for these two PNe are plotted in Fig. 1 . 

Fig. 1 shows that in all nebulae and for all ions where ionic 
abundances from both ORLs and CELs have been derived, the val- 
ues from the optical recombination lines are consistently higher 
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than those derived from the colhsionally excited lines (UV, opti- 
cal or IR). This includes the Magellanic Cloud nebulae SMC N87, 
LMC N66 and LMC N141 ; no such comparisons had been reported 
so far regarding extragalactic PNe. In Table 9 we present a listing 
of the ionic abundance discrepancy factors, ADF(X'+), defined as 
the ratio of the ORL ionic abundance to the UV, optical (OPT) or 
IR GEL ionic abundance. For comparison, data for the previously 
studied planetary nebulae NGC 7009 (LSBC), NGC 6153 (Liu et al. 
2000), the galactic bulge nebulae M 1-42, M 2-36 (Liu et al. 2001b), 
as well as for NGC 6644 and NGC 6572 (from unpuWished ESO 
1.52 m observations) are also included. 

Possible exceptions to the general pattern of large ADF(X*+) 
values include C'''+/H+ in NGC 3918, where abundances derived 
from the C III AA4187, 4650 ORLs and the C IV A1550 reso- 
nant doublet almost coincide; for the same nebula the N*+/H+ 
abundance obtained from the N IV A4606 ORL is 10 per cent less 
that the value derived from the N V A1240 CEL; also N^+/H+ in 
NGC 6302, where identical values are derived from the N ill A4379 
ORL and the N IV] A1486 intercombination line. This could be due 
to an underestimation on our part of the actual electron temperature 
pertaining to highly ionized species in these two high-excitation 
nebulae; adopting a higher Tc in order to derive C'^"'"/H+, N"^+/H+ 
and N*"'"/H"'" would result in lower CEL abundances and produce 
a discrepancy with the ORL values as well. Alternatively, it could 
mean that at least in some nebulae the mechanisms that are respon- 
sible for the abundance discrepancies do not affect species of dif- 
ferent ionization degree in the same way. 

The current work has revealed two further nebulae, NGC 2022 
and LMC N66, which exhibit extreme ORL/CEL ADF's: factors of 
16 and 11, respectively, for O^^fH^. These nebulae are added to a 
rare class of PNe, along with NGC 6153, Ml-42, and Hf 2-2, that 
exhibit ORL/CEL ADF's of ten or more. We also reveal another ob- 
ject, the Type-1 bipolar nebula NGC 2440, whose ORL/CEL abun- 
dance discrepancy factor of 5 is similar to those of NGC 7009 and 
M2-36 (see Table 9). 

Taking into account that in Paper I we presented both 'low- 
density' (adopting [O III] 88 /ini/52 /im densities) and 'high- 
density' values (adopting densities from [Ar iv], [CI ill]), we note 
the following regarding the N^"'"/H+ and O^^/H^ fractions derived 
from the far-IR FS lines that were plotted in Fig. 1: for NGC 3132 
and IC 4406, which have a low mean density, the 'low-density' val- 
ues are plotted for both N^+ and O^"'", while for all remaining PNe 
(including NGC 2022 and NGC 6818), which are of higher mean 
densities, the 'high-density' values are plotted. From our Paper I 
discussion of the results of Rubin ( 1 989) it follows that the low crit- 
ical densities of the far-lR N^"'" and 0'^+ lines lead to their emission 
being biased towards lower-than-average density nebular regions, 
while this is not true for the UV and optical lines, whose emissivity 
ratios relative to H/3 are less sensitive to variations in nebular den- 
sity, owing to the significantly higher critical densities of the lines 
in question. However, it is obvious from Fig. 1 and Table 9 that 
the N^^ and 0"^+ abundances deduced from the IR FS lines agree 
very well with the values obtained from the N ill] A1750 and [O ill] 
AA4959, 5007 lines. This shows that even though density variations 
are present in the nebulae, no significant bias in the inferred CEL 
abundances is present, since very satisfactory agreement amongst 
the derived values is achieved (see also the relevant discussion in 
Section 5.1). Alternatively, this could mean that whatever other bias 
there is, it affects all UV, optical and IR CEL abundances in a rather 
similar manner. 

In Sections 5.1 and 5.2 we will look further into the arguments 



for or against the existence of density or temperature fluctuations 
in our present PNe sample. 

3.3 C/O and N/O elemental ratios 

Tables 10 and 11 respectively present the (C^+/0^+, C/O) and 
(j<[2+/q2+ abundance ratios derived from pure ORL and pure 
CEL line ratios. In both tables the last column may contain two val- 
ues: the one before the vertical dash is computed with the total oxy- 
gen CEL abundance found if we use the O^+(A4959+A5007)/H+ 
ionic ratio in the icf method; the value after the dash adopts the 
oxygen abundance found using the 0^+(A1663)/H+ ratio. In the 
case of the CEL N/O ratio the value after the dash is for ni- 
trogen and oxygen CEL abundances computed by adopting the 
N^+(A1751)/H+ and 0^+(A1663)/H+ ionic ratios, respectively, in 
the icf method. 

The results presented in Tables 10 and 1 1 address an issue that 
we have noted in earlier papers: although CELs and ORLs can often 
yield very different abundances for the same ion relative to hydro- 
gen, ORL-based and CEL-based abundance ratios for two different 
ions are usually quite comparable. This is because ORL/CEL abun- 
dance discrepancy factors are not very different from one ion to the 
next for a given nebula (Table 9). It appears that fairly reliable abun- 
dance ratios for ions (e.g. C^+/0^+) and for elements (e.g. C/O) 
may be obtained, provided that both abundances are based on the 
same type of line (e.g. ORL/ORL or CEL/CEL). It should be noted 
however that any residual differences between such ORL/ORL and 
CEL/CEL elemental ratios may be real and could yield clues about 
the relative nucleosynthetic histories of the ORL and CEL emitting 
media. We urge further study of this issue. In Fig. 2 we show the 
CEL C^+/0^+ and N^+/0^+ ratios compared with the correspond- 
ing ratios derived from ORLs. 

The criterion of Kingsburgh & Barlow (1994) states that Type 
I PNe are those objects that have experienced envelope-burning 
conversion to nitrogen of dredged-up primary carbon, i.e. PNe with 
a N/O ratio greater than the (C + N)/0 ratio of H II regions in the 
host galaxy, the latter being 0.8 for our own Galaxy. According to 
this criterion three objects in our sample, NGC 2440, NGC 6302 
and My Cn 18, all bipolar, qualify as Type 1 PNe by virtue of their 
ORL N/O ratios; but only the first two qualify based on their CEL 
N/O ratios. 



4 CORRELATIONS 

4.1 ORL/CEL ADF's versus Te discrepancies 

The 0^+ ORL/CEL abundance discrepancy factor is known to be 
strongly correlated with the difference between the [O III] forbid- 
den line and H I Balmer jump electron temperatures (Liu et al. 
2001b). In Fig. 3 (top) we plot the ratio of C'^+/H+ ionic abun- 
dances derived from the C ii A4267 ORL to those from the C ill] 
A 1908 CEL versus the difference between the [O iii] nebular to au- 
roral forbidden line temperature and the Balmer jump temperature, 
for a sample of 17 planetary nebulae. 
We see that a strong, positive Unear correlation exists between 

ADF(C2+/H+) = l0g(C2+/H+)oRL - log(C2+/H+)Al908CBL, 

and 

AT = Ta([0 III]) - re(BJ). 

A linear fit to the 17 PNe plotted in Fig. 3 (top) yields. 
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Figure 1. Comparison of ionic abundances derived from optical recombination lines, and from UV, optical and IR coUisionally excited lines. In these and 

subsequent graphs the plotted values are for the entire nebulae originating from ground-based scanning optical spectroscopy and wide aperture, lUE, ISO 
and IRAS observations, except for NGC 6302, NGC 6818, IC4406 and MyCn 18, where abundances from ORLs and optical CELs have been derived from 
fixed-sUt spectra. 
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Figure 1. — continued. 
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Table 9. Ionic ORL/CEL abundance discrepancy factors. 



PN ADF(C2+) ADF(C3+) ADF(N^+) ADF(N3+) ADF(N-4+) ADF(0^+) ADF(03+) ADF(Ne2+) 

ORL/UV ORLAJV ORLAJV ORLAR ORLAJV ORL/UV ORL/OPT ORL/UV ORL/IR ORLAJV ORL/OPT 



NGC2022 


10.(9.1)" 


4.7 


34. 


11. 


15. 


* 


16. 


23. 


9.0 


1.3 


33. 


NGC2440 


4.0(2.5) 


4.8 


6.7 




5.3 


14. 


5.4 


8.9 


* 


* 


2.9 


NGC3132 


4.4 


* 


3.4 


3.9 


* 


* 


2.4 


3.5 


2.8 


* 


* 


NGC 3242 


5.2(5.5) 




7.8 


4.1 






2.2 


5.7 


2.0 


* 


2.5 


NGC3918 


2.1(3.4) 


1.1 


2.3 


1.6 


9.5 


0.9 


1.8 


2.3 


1.4 


* 


8.7 


NGC 5315 


3.4 




1.8 


1.2 






2.0 


1.4 


1.7 




1.3 


NGC 5882 


2.9 


* 


1.5 


2.2 


* 


* 


2.1 


2.2 


2.1 


* 


2.3 


NOP fi'^O? 


6 5 n 

U.J ^J. J J 


* 


5.5 


1.8 


1.0 


2.4 


3.6 


3.5 


3.3 


* 


36. 


NGC 68 18 


2.3(1.3) 


6.2 


12. 


3.7 


14. 


2.3 


2.9 


4.9 


2.1 


1.3 


13. 


IC4191 


* 


* 


* 


3.7 


* 


* 


2.4 


* 


1.5 


* 


3.0 


IC4406 


3.7 


4.9 


1.7 


2.1 


* 


* 


1.9 


1.8 


1.8 


* 


* 


MyCnlS 




* 


* 


* 


* 


* 


1.8 


* 


* 


* 


* 


SMCN87 


2.1 


3.5 






* 


* 


2.8 


* 


* 


* 


* 


LMCN66 


* 


* 


* 


* 


* 


* 


11. 


19. 


* 


* 


* 


LMCN141 


4.6 


2.9 


* 


* 


* 


* 


2.6 


* 


* 


* 


* 


NGC 4361 


15. 


4.0 








* 












NGC 6153 


9.0(9.7) 




3.3 


9.0 






9.2 




7.3 




11. 


NGC 6572 


2.1(1.1) 


* 


* 




* 


* 


1.5 


* 


3.9 


* 


* 


NGC 6644 


2.8 (2.8) 


* 


* 


* 


* 


* 


1.9 


* 


* 


* 


* 


NGC 7009 


4.1(5.3) 


3.2 


3.2 


7.1 


* 


* 


5.0 


* 


5.9 


* 


* 


Ml-42 


23. 


* 


* 


8.2 


* 


* 


22. 


* 


11. 


* 


17. 


M2-36 


4.8 


* 


* 


5.6 


* 


* 


6.9 


* 


5.6 


* 


7.7 



" Parenthesized values are from Rola & Stasiliska (1994). 



Table 10. Comparison of C/0 abundance ratios from CELs and ORLs. 



PN 










ORLC/O 


CEL C/0 




A4267/0 11 


A1908/A1663 


A1908/A4959+A5007 


A 1908/52+88 /jm 






NGC 2022 


0.68 


1.52 


1.08(1.21)'' 


0.60 


0.82 


0.47 


NGC 2440 


0.85 


1.93(1.23)" 


1.16(0.94) 




0.99 


0.96 


NGC 3132 


0.81 


0.64 (0.63) 


0.44(1.68) 


0.51 


0.81 


0.48 


NGC 3242 


0.98 


1.06 (0.66) 


0.42 (0.62) 


0.37 


1.01 


0.4211.05== 


NGC 39 18 


0.94 


1.03(1.15) 


0.83 (0.90) 


0.61 


0.60 


0.60 


NGC 53 15 


0.77 


0.32(0.31) 


0.45 (0.27) 


0.39 


0.58 


0.34 


NGC 5882 


0.39 


0.30 


027 


0.27 


0.44 


0.27 


NGC 6302 


0.61 


0.33 (0.38) 


0.34(0.39) 


0.31 


0.75 


0.31 


NGC 6818 


0.64 


1.35(1.45) 


0.81 (1.28) 


0.60 


0.70 


0.49|0.64 


IC4191 


0.40 


* 


* 


* 


0.45 


* 


IC 4406 


1.09 


0.55 


0.58 


0.54 


1.25 


0.62 


MyCnl8 


0.59 


* 


* 


* 


0.59 


* 


SMCN87 


2.28 


* 


3.04 


* 


2.52 


3.49 


LMCN66 


* 


0.13 


0.08 


* 


* 


0.10 


LMCN141 


1.38 


* 


0.77 


* 


1.48 


1.02 


Solar 


* 


* 


* 


* 


* 


0.50 <^ 



" Parenthesized values are from Rola & Stasihska (1994); 
Parenthesized values are from Rola & Stasinska (1994); 

See Section 3.3 for explanation when two values are Usted for the CEL C/0 ratio; 

Solar oxygen and carbon abundances taken from Allende Prieto et al. (2001, 2002), respectively. 



q2+ 

ADF(— — ) = (0.419 ± 0.051) + (14.6 ± 1.8) x 10"^ AJ, (1) 

with a linear correlation coefficient of 0.91. 

The data include those for seven PNe from other stud- 
ies: NGC 7009 (LSBC), NGC 4361 (Liu 1998; who quotes 
C'^^(A1908)/H+ values from Torres-Peimbert, Peimbert & Pena 
1990), NGC 6153 (Liu et al. 2000), M 1-42 and M2-36 (Liu et al. 
2001b), NGC 6644 (our unpublished observations), and NGC 2392 



[ADF(C^+) from Barker 1991; Tc(BJ) from Liu & Danziger 
(1993b)]. The electron temperatures for the PNe in the current sam- 
ple are taken from Table 7 of Paper 1. 

Fig. 3 (bottom) shows that a very similar linear corre- 
lation exists between ADF(0^+/H+) = log(0^+/H+)oRL - 
log(O^+/H+)A4969+B007 CEL and AT too (cf. Liu et al. 2001b for a 
similar fit to a subset of the nebulae discussed here). A linear fit to 
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Table 11. Comparison of N/0 abundance ratios from CELs and ORLs. 



PN 


N2+/02+ 


N2+/02+ 


N2+/02+ 


N2+/02+ 


ORLN/0 


CEL N/0 




N ll/O 11 


A1750/A1663 


A1750/A4959+A5007 


57 (im/52+88 (im 






NGC2022 


0.34 


0.23 


0.16 


0.27 


0.38 


0.13 


NGC2440 


0.91 


1.21 


0.73 




0.99 


0.59 1 1.09" 


NGC3132 


0.42 


0.42 


0.29 


0.30 


0.51 


0.36|0.43 


NGC 3242 


0.27 


0.20 


0.08 


0.13 


0.29 


0.10|0.26 




u.zo 


U.Z / 


U.ZZ 


u.zz 


fl 1 


U.Zj 


NGC 5315 


0.40 


0.32 


0.45 


0.56 


0.47 


0.54|0.42 


NGC 5882 


0.23 


0.25 


0.23 


0.22 


0.29 


0.33 


NGC 6302 


1.72 


1.10 


1.12 


3.16 


1.53 


1.34 


NGC 6818 


0.58 


0.22 


0.14 


0.34 


0.58 


0.16|0.21 


IC4191 


0.39 


* 


* 


0.16 


0.56 


0.06|0.17 


IC4406 


0.27 


0.28 


0.29 


0.23 


0.34 


0.36|0.34 


My Cnl8 


3.16 


* 


* 


* 


>2.0 


0.39 


SMC N87 












0.33 


LMCN66 




0.55 


0.34 


* 


* 


0.41 


LMCN141 


* 


* 


* 


* 


* 


0.45 


Solar'' 


* 


* 


* 


* 


* 


0.19 



° See Section 3.3 for explanation when two values are listed for the CEL N/0 ratio; 

Solar oxygen abundance taken from Allende Prieto et al. (2001); solar nitrogen abimdance from Grevesse, Noels and Sauval (1996). 



the sixteen planetary nebulae plotted in the bottom panel of Fig. 3, 
yields, 

ADF(-— ) = (0.224 ± 0.053) + (21.6 ± 2.4) x 10"^ AT,(2) 

with a linear correlation coefficient of 0.92. 

The correlations of the abundance discrepancy factors for dou- 
bly ionized C and O with the difference between the [O ill] and 
B aimer jump nebular temperatures are significant, since they pro- 
vide a strong observational indication that the nebular thermal 
structure is intimately tied to the problem of discordant — ORL ver- 
sus CEL — abundance determinations. 



4.2 Abundance discrepancy factors versus PN intrinsic 

surface brightnesses 

Nebular surface brightoess can be considered as an evolutionary 
parameter, since for an expanding nebula it decreases as the neb- 
ula ages and the mean density drops. Gamett & Dinerstein (2001a) 
found that for a sample of about a dozen PNe the magnitude of the 
ADF(0^^/H'^) discrepancy factor was inversely correlated with 
the mean nebular Balmer line surface brightoess, 'suggesting that 
the abundance problem is a function of PN evolution'. In Fig. 4 
we have plotted ADF(02+/H+) and ADF(C2+/H+) values, taken 
from the same sources as used for Fig. 3 (described in Section 4. 1), 
against mean nebular H/3 surface brightnesses, for samples of 21 
and 20 PNe, respectively. Here the surface brightness, S(H/3), is de- 
fined as the flux received per square arcsec of the nebula corrected 
for interstellar extinction. For all galactic PNe we use nebular an- 
gular radii, H/3 integrated fluxes and logarithmic extinction coeffi- 
cients, ciUPY""^, from Cahn, Kaler & Stanghellini (1992; CKS92); 
for the Cloud PNe integrated fluxes are from Meatheringham et al. 
(1988); for LMC N66 we adopt an angular radius of 1.5 arcsec 
(Dopita et al. 1993), for LMC N141 we adopt 0.30 arcsec (Shaw et 
al. 2001), while for SMC N87 we adopt 0.23 arcsec (StanghelUni 
et al. 2003). The ORL/CEL abundance ratios for NGC 6572 were 
derived from our unpublished ESO 1.52-m observations. 



From Fig. 4, we confirm that the abundance discrepancy fac- 
tors are inversely correlated with the nebular surface brightoess. A 
linear fit to the 21 PNe plotted in Fig. 4(top) yields, 

ADF(-^) = (-3.52 ±1.41)- (0.327 ± 0.113) X log ^(H/?), (3) 

with a linear correlation coefficient of —0.55; while a Unear fit to 
the 20 PNe plotted in Fig. 4(bottom) yields, 

ADF(-— ) = (-2.88 ±1.00) -(0.285 ±0.080) x logS'(H/?), (4) 

with a linear correlation coefficient of —0.64, i.e. the fit and cor- 
relation coefficient is not very strong for ADF(0^^/H"'") (Fig. 4, 
upper), being quite better for ADF(C^"'"/H+) (Fig. 4, lower). 



4.3 Abundance discrepancy factors versus PN absolute radii 

Fig. 5 plots the ADF(0^+/H+) and ADF(C^+/H+) abundance 
discrepancy factors versus nebular absolute radii. For all galactic 
PNe we have used the absolute radii quoted by CKS92, except for 
NGC 3132 where the radius was derived assuming a distance to the 
nebula of 600 pc (Sahu & Desai 1986) and an angular radius of 
22.5 arcsec (CKS92). For the Cloud PNe we adopted the follow- 
ing: for LMC N66 a radius of 0.324 pc (Dopita et al. 1993); for 
LMC N141 a radius of 0.074 pc, adopting a nebular angular radius 
of 0.30 arcsec (Shaw et al. 2001) and assuming the same distance to 
the LMC (50.6 kpc) as adopted by Dopita et al. (1996); finally for 
SMC N87 a radius of 0.067 pc, based on the (photometric) optical 
HST diametre of 0.45 arcsec (Stanghellini et al. 2003) and assum- 
ing a distance to the SMC of 60 kpc. 

We see that a positive linear correlation exists between 
ADF(C^+/H+) and the absolute nebular radius for 20 PNe 
[Fig. 5(bottom)], which can be fitted by, 

ri2+ 

ADF(-— ) = (0.242 ± 0.121) -|- (5.51 ± 1.35) x R, (5) 
with a linear correlation coefficient of 0.69. 
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Figure 1. — continued. 
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Figure 2. Comparison of the C^+/0^+ (top) and N^+/0^+ (bottom) ionic 
abundance ratios derived from optical recombination lines, and from UV, 
optical and far-IR coUisionally excited lines. Data points with identical x- 
axis values share the same label. The dashed lines have a slope of unity. 



No such correlation is obvious regarding ADF(0^^/H"'") 
[Fig. 5(top)]. 



4.4 Temperature differences versus PN radii and surface 
brightnesses 

So far we saw that abundance discrepancy factors are: i) positively 
correlated with the difference between the [O iii] forbidden line 
and BJ temperatures, ii) weakly correlated with decreasing intrin- 
sic nebular surface brightness, and iii) positively correlated with the 
absolute PN radii [for the case of ADF(C^+/H+) mostly]. There- 
fore it is not surprising to find that the aforementioned temperature 
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Figure 3. (top): ORL/CEL ionic abundance discrepancy factors for C'^+ 
plotted against tlie difference between the [O III] and BJ electron temper- 
atures; (bottom): the same for the 0^+ ratios; the soUd lines are the linear 
fits of Eqs. 1 and 2; see text for details. 



difference is also inversely correlated with 5(11/3) and positively 
correlated with the absolute nebular radii (Fig. 6). 

A linear fit to the 25 nebulae plotted in Fig. 6 (left) yields, 

Ar= (-17.6 ± 6.1) X 10^ - (1517 ± 481) x log ^(H/?), (6) 

with a linear correlation coefficient of —0.55. 

The relation between AT and the nebular radii for the same 
25 objects [Fig. 6 (right)] can be fitted by, 

AT= (-853 ± 744) + (28.1 ± 7.6) x 10^ R, (7) 

which has a linear correlation coefficient of 0.61. 

In Fig. 6, apart from the 1 1 nebulae whose B J temperatures 
were measured in the context of this study, we also included data 
for 6 objects published previously (see Section 4.1 for details on 
these), along with data for a further 8 nebulae whose BJ and [O iii] 
temperatures were presented by Liu & Danziger (1993b). 



Figure 4. 0^+/H+ (top) and C'^+/H+ (bottom) abundance discrepancy 
factors, plotted against nebular H/3 surface brightness. The soUd lines are 
the hnear fits discussed in the text. 



5 THE EFFECTS OF NEBULAR DENSITY AND 
TEMPERATURE VARIATIONS 

5.1 Density inhomogeneities 

In Paper I nebular electron densities derived from various diagnos- 
tic ratios were presented. It was shown that the low critical den- 
sity [O III] 52- and 88-/im lines yield electron densities that are, 
on average, a factor of 6 lower than those derived from the optical 
[Ar IV] and [CI III] doublets, which have much higher critical den- 
sities. Given that the [O III], [Ar IV] and [CI III] lines arise from 
regions of similar ionization degree,'^ this result points to the pres- 
ence of strong density variations within the nebulae. The effects 
of such variations on the G^+ZH"'" and N^+ZH"*" abundance ratios 
derived from the far-IR lines were discussed in Paper I in the con- 

^ The ionization potential of 0+ (= 35.1 eV), falls between those of C1+ 
(= 23.8 eV) and Ar2+ (= 40.7 eV). 
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Figure 6. The difference between the nebular electron temperatures derived from the [O iii] nebular to auroral forbidden line ratio, Te([0 iii]), and from the 
nebular Balmer continuum discontinuity, Te(BJ), is plotted against the intrinsic nebular H/3 surface brightness (Ze/t) and absolute radius (right). Open circles 
data are from LD93b. The solid lines are the linear fits discussed in the text. 



text of Rubin's (1989) examination of the effects of varying A^c 
on nebular abundances; the general conclusion from Rubin's work 
is that when there are variable density conditions, the intensity ra- 
tios of forbidden lines with low critical densities relative to H/3 
(or any recombination line with a volume emissivity oc Ne^), will 
generally underestimate the true ionic abundances. We saw that in 
accord with these theoretical predictions, abundances from far-IR 
lines are, in general, underestimated relative to UV and optical val- 
ues (that are subject to smaller bias due to their higher A''cr), when 
the 'low' [O III] 52/im/88/im densities are adopted for the calcula- 
tions. On the contrary, we saw that when higher A'^o's were adopted 
(intermediate or equal to the [Ar iv] and [CI ill] densities), the in- 
ferred abundances from the far-IR lines become compatible with 
those derived from the UV and optical indicators (see Fig. 1 of this 
paper and Fig. 3 of Paper I); with the exception of the high critical 
density [Ne III] 15.5-/im line which is not affected by such modest 
density variations. 

In all cases, however, and even after accounting for mod- 
est density variations, the CEL abundances derived from UV, op- 
tical or IR lines remain consistently lower than the ORL values 
(Fig. 1). Viegas & Clegg (1994) showed that dense clumps with 
A^o > 10^' cm^'^ can have a substantial effect on the derived [O III] 
forbidden-line temperature, via coUisional suppression of the neb- 
ular AA4959, 5007 lines; in this way the observationally derived 
Tc([0 III]) could be significantly overestimated and CEL abun- 
dances accordingly underestimated. In their analysis of the extreme 
nebula NGC 6153, Liu et al. (2000) discussed the case for the ex- 
istence of dense condensations with a very small filling factor as a 
potential explanation for the large ORL/CEL abundance discrepan- 
cies. They concluded that such dense clumps would have a different 
effect on different CEL lines. In particular, they found that in the 
increased density environment of the condensations, most IR and 
optical CELs would be collisionally suppressed, accounting for the 



abundance ratios, thus covering part of the distance from the ORL 
values; on the other hand however, the abundances derived from 
the high critical density UV CELs, using the downward corrected 
electron temperatures that result from allowing for the coUisional 
suppression of the AA4959, 5007 lines, would become so enhanced 
as to exceed the ORL abundance results. The conclusion in this 
case is that the existence of dense condensations is not in itself a 
sufficient solution to the problem, since it requires the discrepan- 
cies to be correlated with the critical densities of the various CELs, 
something not indicated by the observations. 

5.2 Temperature fluctuations 

Peimbert (1967) first proposed that in the presence of tempera- 
ture fluctuations within nebulae, adoption of the [O III] (A4659 
+ A5007)/A4363 line ratio as a standard thermometer would re- 
sult in the underestimation of elemental abundances as derived 
from CELs, since that line ratio would be biased towards high- 
temperature nebular regions. Generally speaking, if variations in 
electron temperature exist along the line of sight through a nebula, 
abundance ratios derived from the ratio of a collisionally excited 
line to a recombination line (e.g. the most commonly used recom- 
bination line, H/3) will be underestimated, while those derived from 
a pure recombination line ratio should be almost insensitive to Te 
and thus largely unaffected by any incorrect assumptions regarding 
the nebular temperature distribution. 

Peimbert (1971) found possible evidence for the existence of 
such temperature fluctuations, by comparing for three PNe the tem- 
perature derived from the ratio of the Balmer continuum jump to the 
intensity of H/3 with that obtained from the [O III] nebular to auro- 
ral CEL ratio; he found that the latter temperatures were higher 
than Te(BJ). This result has been supported by further observa- 
tions of more planetary nebulae (e.g. Liu & Danziger 1993b). The 
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Figure 5. 0^+/H+ (top) and C-^+/H+ (bottom) abundance discrepancy 
factors, plotted against nebular radius. The solid line in the lower plot is the 
linear fit discussed in the text. 



Balmer jump temperatures that we presented for 12 PNe in Table 7 
of Paper I were also systematically lower than the corresponding 
re([0 III]) values presented there. 

The concept of temperature fluctuations in a chemically ho- 
mogeneous nebula has been invoked several times as a promising 
cause of the ORL/CEL abundance discrepancy problem, pertaining 
to both discrepant C^+(A4267)/H+ and 0^+/ll+ ORL abundance 
ratios. For instance, Peimbert et al. (1993) proposed that spatial 
temperature fluctuations with an amplitude of ~ 20 per cent around 
a mean value would reconcile the factor of 1.55 discrepancy be- 
tween the ORL and CEL abundances in the planetary neb- 
ula NGC6572. Peimbert, Luridiana & Torres-Peimbert (1995) de- 
rived electron temperatures for a sample of nebulae from the C ill] 
A1908/C II A4267 ratio and showed that they are generally lower 
than Tc([0 III]). On account of this they argued that in order to ex- 
plain the discrepancy between the C^+(A4267)/O^+(A5007) and 
C^+(A1908)/O^+(A5007) values in planetary nebulae published 
by Rola & Stasinska (1994), large temperature fluctuations were 



needed; they also recommended the adoption of re(C ) rather 
than To([0 III]) when deriving ionic abundances from CELs and 
concluded that due to their insensitivity to temperature variations, 
abundances derived from pure recombination line ratios are more 
reliable. 

However, the validity of the C ill] A1908/C ii A4267 ratio (or 
even that of the C iv A1549/C iii A4650 ratio) as a reliable nebular 
thermometer is ambiguous, since there is evidence that, at for PNe 
which exhibit non-uniform, inhomogeneous abundances, the C II 
A4267 ORL and C III] A 1908 CEL do not originate in the same 
gaseous component. Harrington & Feibelman (1984) argued that 
such is the case for the hydrogen-deficient material which exists at 
the centre of the planetary nebula A 30 in the form of condensa- 
tions. Subsequent, detailed photoionization modelling of HST and 
ground-based spectra of several of the knots in A 30 (Borkowski et 
al. 1993, Ercolano et al. 2003) has confirmed that the ORLs and 
CELs originate from different regions of the knots. 

We have investigated our results for evidence that temperature 
fluctuations may contribute to the observed abundance discrepan- 
cies. The major premise of this scenario is that ionic abundances 
derived from the intensity ratio of two lines with very different 
temperature dependencies — e.g. 7([0 ill] A1663)//(H/?) — will be 
severely biased if an incorrect Tc is adopted in calculating line 
emissivities. Since, on average, re(BJ) < re([0 ill]), and since the 
standard practice is to use the [O iii] optical Une ratio temperature 
in order to calculate both the emissivities of CELs and of H/3, the 
possibility of systematically underestimating CEL abundances is 
there; should then ORL abundances that are insensitive to tempera- 
ture bias be trusted more and if so why do they point towards heavy 
element overabundances in most nebulae? One way to try to answer 
this question is to examine the effects of changing the adopted Tc 
on ionic abundances derived from CELs that have different temper- 
ature dependencies. 

In the current study the standard electron temperature from 
the [O III] (A4659 -i- A5007)/A4363 line ratio was adopted in or- 
der to derive all ionic abundances of doubly ionized species. If we 
suppose however that the Te pertinent to the nebular zones where 
doubly ionized species exist is lower than that, e.g. that it may 
be as low as the Hi BJ temperature, then the ionic abundances 
derived from CELs would all be underestimated by our adoption 
of Te = rc([0 III]), but not all by the same factor. Abundances 
derived from high excitation energy transitions, such as O ill] 
A1663 (Sex ^ 86 kK), N iii] A1751 (Bex ^ 82 kK), or C ill] A1908 
(Bex — 75 kK), would be underestimated more relative to those de- 
rived from transitions that have a lower temperature sensitivity, 
such as the [O III] AA4959, 5007 nebular lines (Ec^~29kK). 
Our results show however (Table 9), that all CEL C^""", N^""" and 
0^+ abundances are systematically lower than those derived from 
ORLs, by rather similar factors. 

For instance, the average C^+ A4267/A1908 ADF for eigh- 
teen nebulae is 5.4, while the 02+(ORL/A4959 -I- A5007) ADF 
is 5.1; also flie 0^+ ORL/AI663 ADF for ten nebulae is 5.7 (Ta- 
ble 9). This apparent uniformity caimot be reconciled with the dif- 
ferent behaviour that one would expect amongst the discrepancy 
ratios in the case of temperature fluctuations of the type envis- 
aged by Peimbert. We would expect the C^+ ORL/CEL ADF's to 
be greater than that for 0^+, since the 75 kK excitation energy 
of the upper levels of the A 1908 doublet exceeds the 62 kK ex- 
citation energy of the upper level of the A4363 transition. How- 
ever, in flie extreme nebula NGC 2022 the ADF for C^+ is 10, but 
that for 0^+ is 16, i.e. the opposite of that expected; in the other 
extreme nebulae NGC 6153 (Liu et al. 2000) and M 1-42 (Liu et 
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al. 2001b) the corresponding factors are almost equal. Similarly, 
we would expect the C^^ ORL/CEL ADF to be greater than that 
for C^"*", due to the higher excitation energy of the C IV A1549 
resonant doublet (Sex~93kK), relative to that of C iii] A1908 
(_Eex — 75 kK); however, this is not the case for the nine PN with 
both A4267/A1908 and ORL/A1549 ADF's documented, 
where on average, these abundance discrepancy factors are 5.3 and 
3.9, respectively. 

Having just seen that in the presence of temperature fluctua- 
tions we would expect ADF(C^^/H^), as defined in the previous 
section, to be affected more by errors in Te than in ADF(0^+/H"''), 
our fits (Eqs. 1 and 2) show that actually it may well be that the 
opposite is true, judging from the slightly steeper ADF(0^^/H+) 
function (Fig. 3). This behavioiu' of the two discrepancy factors 
conflicts with that expected for temperature fluctuations of the stan- 
dard type. 

Probably, however, the most serious obstacle to the classical 
temperature fluctuation hypothesis lies with the IR abundance re- 
sults. Inspection of Fig. 1 and Table 9 shows that the ORL/CEL 
ADF's for CEL abundances obtained from the far-lR [O III] and 
[N III] lines are in Une with those from the optical and UV CELs, 
i.e. there is no correlation whatsoever with the very low i5ex 
(< 1000K) of IR transitions. As we mentioned in Section 3.2, 
very satisfactory agreement is reached amongst UV, optical and IR 
CEL 0^+/H+ ADF's, while the same is also true for N^+/H+ and 
Ne^+ZH"*" ; this would not be the case in the presence of strong tem- 
perature fluctuations in a chemically homogeneous nebula. 

Ruiz et al. (2003) and Peimbert et al. (2004) have presented 
deep optical echelle spectra of NGC5307 and NGC5315, respec- 
tively, deriving 0^+ ADF's of 1.9 and 1.7 for these two PNe (we 
obtained an optical ADF of 2.0 for NGC5315; see Table 9). 
They argued that these ADF's could be explained by classical 
temperature fluctuations, corresponding to values of 0.051 and 
0.056, respectively. NGC 5307 has no measurements available for 
its far-IR FS lines, but ISO spectra exist for NGC 5315. Peimbert 
et al. (2004) chose not to use the the 52- and 88-^tm [O III] ISO 
LWS line fluxes presented by Liu et al. (2001a), on the grounds that 
the electron density of 2290 cm~^ yielded by the intensity ratio of 
these two Unes lies between the respective critical densities of 3800 
and 1800 cm"^. However in Paper I we found that for NGC 5315 
(and for all the other PNe in our sample with electron densities ex- 
ceeding one or both of these critical densities), the adoption of an 
electron density given by the mean of the densities from the [CI ill] 
and [Ar iv] optical doublet ratios gives an 0^+ abundance from the 
far-IR fine-structure lines that agrees with the value obtained from 
the optical forbidden lines. 

With regard to this issue, two PNe in our current sample, 
NGC 3132 [ADF(0^+) = 2.4] and IC4406 [ADF(02+) = 1.9], are 
of key importance, since their far-IR [O III] doublet ratio densities 
of 355 and 540 cm^"^ (Paper I) are well below the critical densi- 
ties of the two far-IR lines (as are the electron densities derived for 
these two nebulae from the various optical doublet ratios; see Ta- 
ble 6 of Paper I). Thus there should be no objection to the derivation 
of O^^ abundances from the far-IR [O III] lines for these nebulae. 
Inspection of Table 9 of Paper I shows that for NGC 3132 the 
abundance derived for the far-IR FS lines, using the far-IR doublet 
ratio density, differs by less than 14% from the O^""" abundance de- 
rived from the optical forbidden lines, while for IC 4406 they differ 
by less than 3%. 

Expressing this another way, the electron temperature of 
9900 K derived for both NGC 3132 and IC4406 from the ra- 
tio of their IR to optical [O ill] doublet fluxes agrees with the 



temperatures of 9530 K and 10000 K derived from their [O III] 
A4363/(A5()07+A4959) ratios (Table 7 of Paper I). If the rather 
typical ADF's of ~2 derived for these two PNe were caused by 
temperature fluctuations in a chemically homogeneous medium, 
then the IR to optical [O III] doublet flux ratio should yield sig- 
nificantly lower temperatures than the much more temperature- 
sensitive A4363/(A5007+A4959) ratio. This is not the case, so we 
conclude that classical temperature fluctuations cannot explain the 
ADF's of ~2 foimd for those two nebulae. 



6 EVIDENCE FOR COLD PLASMA FROM ORL 
ELECTRON TEMPERATURES 

6.1 Electron temperatures from He I ORL ratios 

Liu (2003) used the weak temperature sensitivity of several He I re- 
combination lines to derive He+ electron temperatures for five PNe 
with particularly high ADF's. He found that both the value of 
Te(He I) and the ratio of Te(He I)/Te([0 ill]) decreased system- 
atically with increasing nebular ADF, with Tc(He I) ranging from 
5380 K for NGC 7009 (ADF = 4.7) down to 2310 K for M 1-42 
(ADF = 22) and as low as 775 K for Hf 2-2 (ADF = 84). This 
provided observational support for the presence of regions of cold 
plasma within these high-ADF nebulae. 

It is of interest to investigate what patterns might be shown 
by the lower- ADF PNe in the current sample, so we have derived 
He I electron temperatures using the curves presented by Liu (2003; 
based on the data of Smits 1996). We used his A6678/A4471 curves 
only, since the A5876/A4471 curves have minima between 6000- 
10000 K, which allow matches to a given line intensity ratio at 
both low and high temperature. Column 5 of Table 12 presents the 
dereddened A6678/A4471 flux ratios (from Paper I) and the derived 
values of re(He i) are listed in column 6. Also presented in Table 12 
are the nebular electron densities (column 2; adopted from Ta- 
ble 6 of Paper I) that were used to interpolate between the density- 
dependent A6678/A4471 curves, as well as the nebular [O ill] and 
hydrogen Balmer-jump electron temperatures (columns 3 and 4), 
which were taken from Table 7 of Paper I. Very accurate relative 
flux calibration over the wavelength range spanning the two He I 
lines is required, since the A6678/A4471 ratio varies by only 55- 
60% between electron temperatures of 1000 K and 16000 K. De- 
spite the stringent requirement on the accuracy of the relative flux 
calibration over such a wide wavelength range, six of the nine PNe 
listed in Table 12 show He I electron temperatures within 2000 K 
of the nebular [O III] electron temperatures, one (NGC 3132) has a 
He I electron temperature which is 3300 K larger than the [O ill] 
electron temperature, and two of nebulae, NGC 6818 and IC 4191, 
yield He I electron temperatures that are respectively 8300 K and 
7300 K lower than the corresponding [O ill] electron temperatures. 



6.2 Electron temperatures from O 1 1 ORL ratios 

In this subsection we present observational evidence that the O II 
ORLs detected from several of our PNe are emitted from ionized 
gas which is at temperatures much lower than those derived from 
either the Hi Balmer discontinuity or the [O ill] optical forbidden 
lines, for which values were derived in Paper I (see also the current 
Table 12). We will make use of the temperature sensitivity of the 
ratio of the intensities of the A4089 and A4649 Unes of O ii, a tech- 
nique that was first used by Wesson, Liu & Barlow (2003), who 
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Table 12. Comparison of CEL and ORL electron temperatures (in K). 



Nebula A^e TeCLO III]) Te(BJ) /(A6678) Te(He I) /(A4649) /(A4089) Te(Oll) ADF(0^+) 

(cm-3) CEL BJ/HU //(A4471) A6678/A4471 On /7(4649) A4089/A4649 ORL/CEL 



NGC 2022 


1500 


15000 


13200 


0.693 


15900 


0.333 


0.607 


<300 


16 


NGC 2440 


6000 


16150 


14000 






0.100 


0.420 


<300 


5.4 


NGC 3132 


600 


9530 


10800 


0.725 


12800 


0.195 






2.4 


NGC 3242 


2000 


11700 


10200 


0.774 


10000 


0.216 


0.307= 


2600 


2.2 


NGC 3918 


5000 


12600 


12300 


0.740 


12000 


0.211 


0.288 


3650 


1.8 


NGC 5315 


10000 


9000 


8600 


0.775 


10000 


0.344 


0.261 


5750 


2.0 


NGC 5882 


4000 


9400 


7800 


0.763 


10700 


0.362 


0.230= 


8700 


2.1 


NGC 6302 


14000 


18400 


16400 


0.671 


15100 


0.100 






3.6 


NGC 6818 


2000 


13300 


12100 


0.841 


5000 


0.182 


0.300= 


2900 


2.4 


^4191" 


10000 


10700 


10500 


0.908 


3100 


0.725 


0.284 


3900 


2.4 


IC4191'' 


10000 


10000 


9200 


0.916 


2900 


0.572 


0.316 


2150 


2.4 


IC4406 


1000 


10000 


9350 


0.797 


8000 


0.284 






1.9 



" Fixed-sUt observation; 
' Scanning-sUt observation; 

= Adopted O II 4089.29 A intensity is after correction for Si IV 4088.85 A, using the measured intensity of Si IV A41 16.10 A and 
adopting a 2:1 flux ratio for the 4s2S-4p2p° Si IV A4089/A4116 doublet components. 



derived very low electron temperatures for two hydrogen-deficient 
knots in the planetary nebula Abell 30. 

The intensities of recombination lines originating from states 
of different valence orbital angular momentum have different de- 
pendences on electron temperature. By comparing the intensity of 
a Une in the O ii 4f-3d transition array with one from the 3p-3s 
array, for example, it is possible to deduce a recombination line 
temperature. There are several potential difficulties in using such 
lines to measure the electron temperature. Firstly, the variation of 
the intensity ratios with temperature is weak, typically a factor of 
two or three when the temperature changes from 300 K to 10000 K, 
meaning that it is important to ensure that the intensities of the 
weak recombination lines are determined as accurately as possi- 
ble. Secondly, the relative intensities may be affected by the distri- 
bution of population in the recombining ion ground state, O^"*" 
in this case. The available O ii recombination coefficients assume 
that the levels of this term are populated according to their statis- 
tical weight (Storey 1994). If the electron density is sufficiently 
low, this may not be the case, with the energetically higher levels 
of the O^""" state being relatively underpopulated; this has been 
observed for the low density H ii regions 30 Doradus, LMC Nl IB 
and SMC N66 (Tsamis et al. 2003a). Equally, a very low electron 
temperature would have the same effect. 

To attempt to circumvent this second problem, we consider 
the ratio of intensities of A4089 from the 4f-3d transition array and 
A4649 from the 3p-3s array. These lines originate from the state of 
highest total angular momentum, J, in each case and should there- 
fore both be mainly populated from the O^^ ^P2 level. The pub- 
lished recombination coefficients for these two lines (LSBC and 
Storey 1994, respectively) are only valid for temperatures Tc > 
5000 K. Since we wish to investigate the very low temperature 
regime we have extended the recombination coefficient calcula- 
tions to temperatures Te > 300 K. Figure 7 shows the correspond- 
ing intensity ratio for temperatures between 300 K and 15000 K. 
The calculations could not easily be taken to even lower temper- 
atures due to numerical problems in the recombination coefficient 
codes. 

Measurements of the intensities of A4089 and A4649 were 



given in Paper I. The A4089 line is isolated and measurement 
of the intensity is relatively straightforward.^ The A4649 line, 
on the other hand is part of a complex blend in which we 
identify five components, two lines of the O II multiplet V 1, 
A4649.13 and A4650.84, and three lines of the C III multiplet V 1, 
AA4647.42, 4650.25, 4651.47. We have made new five-component 
fits to this blend for a subset of the PNe, using the known wave- 
length separations between the lines, assuming that all five lines 
have the same width as the C II A4267 line (the strongest heavy 
element ORL in the spectra), and using the expected intensity ratio 
of 5:3: 1 for the three C ill Unes. The resulting A4089 to A4649 Une 
ratios and derived electron temperatures are given in Table 12 for 
the six PNe for which both lines were detected, while the derived 
relative intensities for A4649.13 and A4650.84 can be found in Ta- 
ble B13. Note that the O ii electron temperatures presented here 
in Table 12 supersede those that were listed, without discussion, in 
Table 7 of Paper I. 

Inspection of Table 12 shows that two of the eight PNe 
with derived A4089/A4649 O II electron temperatures show val- 
ues below our 300 K theoretical calculation limit. The fixed-slit 
and scanning-slit observations of IC 4191 yield Te's of 3900 K 
and 2150 K, respectively, while NGC 5882 has the highest mea- 
sured rc(0 II) of 8700 K. It can be shown that these values are 
reduced /wrt/ier when one takes into account the small contribution 
to the O II AA4089, 4649 intensities from the normal nebular gas at 
Te([0 III]) temperatures. Once this contribution is subtracted from 
the observed O II intensities, the revised O II ORL ratios (not shown 
in Table 12) indicate electron temperatures of less than 300 K for 
six PNe, while for NGC 53 15 and NGC 5882 the resulting temper- 
atures are 4350, and 7190 K respectively. These findings point to- 



^ Apart from the cases of NGC 3242, 5882 and 6818, where the detection 
of weak emission in the 4116.1 A component of the Si IV 4s^S^p^P° 
AA4089,4116 doublet (with relative intensities of 0.022, 0.0308 and 0.022, 
respectively, on a scale where H/3 = 100) led us to correct the measured 
4089 A line flux for the presence of Si IV 4088.85 A, by subtracting twice 
the measured intensity of the Si IV 4116.10 A line 
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Figure 7. The theoretical intensity ratio for O II A4089/A4649 as a function 
of electron temperature. See Section 6.2 for details. 

wards the presence of ultra-cold plasma regions in a large fraction 
of our sample PNe. 



7 DISCUSSION 

In Sections 3 and 5.2 it was demonstrated that there is no depen- 
dence of the magnitude of the nebular ORL/CEL abundance dis- 
crepancy factors upon the excitation energy of the UV, optical or 
IR CEL transition used (see Table 9), indicating that classical (i.e. 
in a chemically homogeneous medium) nebular temperature fluctu- 
ations are not the cause of the observed abundance discrepancies. 
This reinforces the same conclusion that was reached in Paper I, 
based there upon the fact that [O III] electron temperatures derived 
from the ratio of the 52 and 88-/im PS lines to the 4959 and 5007- 
A forbidden lines were greater than or comparable to those derived 
from the ratio of the higher excitation energy 4363-A transition to 
the 4959 and 5007 -A lines — if temperature fluctuations in the am- 
bient nebular material were the cause of the ORL/CEL abundance 
discrepancies then the IR line-based ratio should yield lower tem- 
peratures. We conclude instead that the main cause of the abun- 
dance discrepancies is enhanced ORL emission from cold ionized 
gas located in hydrogen-deficient clumps inside the nebulae, as first 
postulated by Liu et al. (2000) for the high-ADF PN NGC 6153. 

When nebular heavy element abundances exceed about five 
times solar, cooling by their coUisionally excited infrared fine- 
structure (IR FS) lines is alone sufficient to balance the photo- 
electric heating from atomic species. Since cooling by the low- 
excitation IR FS lines saturates above a few thousand K, nebular 
electron temperatures in such high-metallicity regions will drop 
to values of this order, the exact value being determined by the 
H to heavy element ratios and the input ionizing spectrum. It is 
therefore physically plausible and self-consistent for the heavy el- 
ement ORLs that yield enhanced abundances relative to those de- 
rived from CELs to also indicate very low electron temperatures for 
the regions from which they emit. 

Liu et al. (2000) in their empirical modelling of NGC 6153 had 
postulated the presence of H-deficient ionized clumps within the 
nebula that would be cool enough to suppress optical forbidden- 
line emission and even some infrared fine-structure line emission 



but which would emit strongly in heavy element recombination 
lines, due to the inverse power-law temperature dependence of 
ORL emission. Pequignot et al. (2002) have constructed photoion- 
ization models of NGC 6153 andM 1-42 incorporating H-deficient 
inclusions and found equilibrium Tc's of ~ lO'' K in the H-deficient 
clumps and ^ 10* K in the ambient gas, with the H-deficient ion- 
ized regions being within a factor of two of pressure equiUbrium 
with the ambient nebular gas. In both cases, the H-deficient model 
components contained only 1% of the total ionized mass, so that 
the overall metallicity of the whole nebula was close to that of the 
'normal' high-temperature component. 

Similar, or lower, mass fractions for the postulated H-deficient 
clumps in the typical PNe studied here, which have lower ADF's 
than the extreme cases discussed above, should ensure that their 
integrated IR FS line emission will not represent a significant per- 
turbation to the integrated IR FS emission from the ambient high- 
temperature nebular material that forms the majority of the neb- 
ular mass, particularly if the H-deficient clump electron densities 
exceed the rather low critical densities of 500 - 3000 cm^^ that 
correspond to the [O III] and [N III] IR FS lines that have been in- 
vestigated in this paper. The strong inverse power-law temperature 
dependence of ORL emission means that material at a temperature 
of 500 K will emit an ORL such as O II A4649 eighteen times more 
strongly than at 10* K. If 0.5% of the nebular mass was located in 
500 K clumps having an electron density ten times that of the ambi- 
ent nebular gas, then the integrated ORL emission from the clumps 
would exceed that from the ambient gas by a factor of nine. 

Hydrogen is not expected to be entirely absent in the postu- 
lated H-deficient clumps. The steep inverse power-law temperature 
dependence of its recombination emission from the clumps could 
be responsible for the strong correlation, discussed in Section 4.1, 
between the C'^^ or O^"'" ADF's and the temperature difference, 
AT, between the [O III] optical forbidden line and H I Balmer 
jump electron temperatures (see Fig. 3). The correlation between 
ionic ADF's and decreasing nebular surface brightness (Fig. 4), 
or increasing absolute nebular radius (Fig. 5), may indicate that 
the density contrast between the H-deficient clumps and the am- 
bient nebular material increases as the nebula evolves, i.e. that the 
clump density decreases less than the ambient nebular density does 
as the nebula expands. Such a situation might arise if the clump ion- 
ized gas originates from photoevaporation of dense neutral cores 
(cometary knots) of the kind found in the Helix Nebula (Meabum 
et al. 1992; O'Dell et al. 2000) and the Eskimo Nebula (NGC 2392; 
O'Dell et al. 2002). Indeed, the location of NGC 2392 in the AT 
vs. SiViP) and AT vs. Radius diagrams (Fig. 6) is consistent with 
this picture and indicates that NGC 2392 is a candidate high-ADF 
nebula. In confirmation of this, optical and ultraviolet large aper- 
ture measurements presented by Barker (1991) for six positions 
in NGC 2392 yielded A4267/A1908 C^+ ADF's ranging from 7 
to more than 24. We thus show NGC 2392 in the ADF(C^+) vs. 
AT diagram (Fig. 3) and the ADV{C'^+) versus Radius diagram 
(Fig. 5) using a mean C^"*" ADF of 18 derived for the positions 
observed by Barker. It would be of interest to obtain deep optical 
spectra of this nebula in order to examine its heavy element ORL 
spectrum in detail. 

Available spatial analyses of long-slit PN spectra show that 
ORL/CEL ADF's peak towards the center of nebulae in the cases of 
NGC 6153 (Liu et al. 2000) and NGC 6720 (Ring Nebula; Gamett 
& Dinerstein 2001b). An examination of the results presented by 
Barker (1991) for NGC 2392 shows that the same trend is present 
in those data too: the inferred ADF(C^"'") increases towards the cen- 
tre of the PN (along the aligned positions 4, 2, and I; cf. Fig. 1 
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of that paper). In NGC6720 especially, the location of peak O II 
ORL emission does not coincide with the positions of the HST- 
resolved dusty cometary knots, which populate the main shell of 
the Ring, but is displaced inwards from that of the peak [O ill] 
emission (Garnctt & Dincrstcin 2001b). We speculate that such an 
effect could be due to the advanced photo-processing of knots in 
the Ring Nebula that were overcome by the main ionization front 
in the past and whose relic material, rich in heavy elements, is now 
immersed in the He^"*" nebular zone, being subjected to the intense 
radiation field of the central star. The question of the possible re- 
lationship between the ORL/CEL abvmdance discrepancy problem 
and the cometary knot complexes observed in many PNe warrants 
further investigation. 

Ruiz et al. (2003) and Peimbert et al. (2004) have ar- 
gued against the presence of H-deficient knots in NGC 5307 and 
NGC5315 on the grovmds that their 34 km s~^ resolution echelle 
spectra did not reveal a difference between the radial velocities or 
line widths of the heavy element ORLs and those of the main neb- 
ular lines, of the type exhibited by the high-velocity H-deficient 
knots in the born-again PNe A30 and A58. We note however that 
the H-deficient knot model that has been invoked to explain the 
ORL/CEL ADF's of typical planetary nebulae makes no specific 
predictions as to whether the knots should exhibit a different kine- 
matic pattern from the bulk of the nebula. Scenarios for the ori- 
gin of such knots in 'normal' PNe include i) the evaporation of 
primitive material (comets, planetesimals) left over from the for- 
mation of the progenitor star (e.g. Liu 2003), which would predict 
ORL C/O and N/O ratios typical of the unprocessed ISM material 
out of which the star formed; or ii) that they originated as incom- 
pletely mixed material brought to the surface by the 3rd dredge- 
up and ejected along with the rest of the AGB progenitor star's 
outer envelope during the PN formation phase. In the latter case, the 
knots' ORL C/O and N/O ratios should show the same nucleosyn- 
thetic signatures as the rest of the ejected envelope. Although the 
postulated H-deficient clumps may not exhibit different kinemat- 
ics from the bulk of the nebular material, the low inferred electron 
temperatures of the clumps and the consequent very low thermal 
broadening means that heavy element ORL's in high-ADF nebulae 
should exhibit much narrower line widths than do the strong forbid- 
den lines. For example, a C ii A4267 line originating from 1000 K 
material should have a FWHM of 2 km s~^ if thermal broaden- 
ing dominates, i.e. if velocity broadening of the line is minimal, 
as when the emission comes from the edge of the nebula, where 
material is moving in the plane of the sky, or comes from well- 
separated approaching and receding velocity components, e.g. near 
the centre of the nebula. Thus observations of suitably chosen neb- 
ular sub-regions at a resolving power of 1.5x10^ may be capable 
of confirming the presence of the postulated cold plasma clumps. 
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APPENDIX A: THE CALCULATION OF ELEMENTAL 
ABUNDANCES FROM RECOMBINATION LINES 

In this Appendix we discuss the derivation of total C, N and O 
abundances for our PN sample from ORL ionic abundances. We 
discuss each PN individually in order to highlight differences in 
the adopted icf scheme. 

NGC3242: the C ii A4267, C ill AA4069, 4187 and A4650 lines 
are detected; we adopt C^+/H+ = 2.03 x lO"**, as an intensity 
weighted mean from the AA4187, 4650 lines only. The unseen C"*" 
is corrected for using the KB94 icf formula: icf{C) = 1 -i- 0+/0^+ 
= 1.01. We calculate the C/H fraction using, 

C/H = ic/(C) X (C^+ -h -I- C*+)/H+. (Al) 

This yields C/H = 8.45x10"'*; this value is only 2 per cent 
more than the one found if C''+/H+ = 1.88 x lO"'', as derived from 
C IV A4658 is not added in. The latter value is an upper limit to the 
C^+/H+ fraction, since A4658 is blended with [Fe ill] A4658.10. 

From an intensity weighted mean of N^"'"/H"'" ratios from 
seven N II lines (cf. Table 4), we find N2+/H+ = 1.71 x 10"". From 
the N III A4379 line, N'^+/H+ = 6.85 x 10"". We use the CEL ratio 
N+/N^+ = 0.0145 (cf. Paper I) to account for the unseen N^; the 
error introduced should be negligible, since N+/N = 0.008. Hence, 
summing three ionic stages, we find N/H = 2.42 x 10"''. 

From our extensive O II recombination-line survey, we derive 
02+/H+ = 6.28 x 10""*, from five O II multiplets (Table 6). An 
ORL 0+/H+ abundance is not available, so the unseen O"*" is cor- 
rected for assuming that the Paper I CEL ratio of 0+/0^+ = 0.0091 
holds for the ORLs too. In view of the minor concentration of O"*" in 
this object (about 1 per cent), the errors introduced should be neg- 
ligible. We then employ the KB94 standard ionization correction 
factor to account for the unseen 0^+, icf{0) = {\+ He^+/He+)^/^ 
= 1.17, and calculate the total O/H abundance, using, 

0/H = ic/(0) x (0+-|-0^+)/H+. (A2) 

This yields, O/H = 7.41 x 10"'*. Alternatively, we can adopt 
03+/H+ = 2.05 X 10"", as derived by LD93a from the O III V8 
multiplet at 3265 A; hence, summing all three ionic stages we find 
O/H = 8.39 X 10"", i.e. 13 per cent larger than yielded by the icf 
method. We adopt this latter value. 

NGC5882: the C ii A4267, C ill AA4069, 4187 and A4650 lines 
are detected; we adopt C^+/H+ = 5.27 x 10"^, as an intensity 
weighted mean from the AA4187, 4650 lines only. No C""*" is ex- 
pected to exist, since only 2 per cent of He is in the form of He^"*" . 
Using Eq. Al with icf{C) = 1.03, we find, C/H = 4.43 x 10"". 

From an intensity weighted mean of five N^+/H+ ratios (Ta- 
ble 4), N^+/H+ = 2.24 X 10"" is found. From the N ill A4379 line, 
N^+/H+ = 5.75 x 10"^ As for NGC 3242, we use the Paper I CEL 



ratio to account for the unseen N+, N+/N2+ = 0.0278 in this case; 
we obtain N/H = 2.88 x 10"". 

From eight O II ORL multiplets we derive 0^+/H+ = 
9.70 X 10"" (Table 6). We account for unseen O"*" as previously, 
using the Paper 1 CEL ratio of 0+/0^+ = 0.029; we have icf{0) = 
1.01, so Eq. A2 yields O/H = 10.08 x 10"". 

NGC 5315: the C ii A4267, C III A4187 lines are detected; no C"+ 
exists since the He^'''/H"'" fraction is negligible. We correct for un- 
seen C+ using the CEL ratio of C+/C^+ = 0.08 (given the small 
concentration of C""", about 7 per cent, no significant uncertainty 
is introduced); hence we sum three ionic stages and find C/H = 
7.29 X 10"". 

From an intensity weighted mean of four N^+/H+ ratios (Ta- 
ble 4), N2+/H+ = 3.43 X 10"" is found; the N III A4379 line is not 
detected. If we assume that N+/N^+ = 0. 163 and N^+/N^+ = 0.565 
ratios, as given by the CELs (Paper I), are also valid for the ORLs, 
we can account for both N"*" and N'^+, hence, N/H = 5.93 x 10"". 
This result may be quite uncertain, since probably the N^+/N^+ 
ratio derived from CELs ls- not equal to the same ratio as derived 
from ORLs (cf. the case of NGC 3918 below). 

From five O II ORL multiplets, 0^+/H+ = 8.57 x 10""; 
although this is a relatively low excitation object, a substantial 
O^"*" concentration has however been derived from the A1401 CEL 
line (Paper 1). If we account for unseen 0+ using the CEL ratio 
0+/0^+ = 0.098, we find O/H = 9.41 x 10"", ignoring 0^+/H+; 
however, if we further assume that 0^"'"/0^"'" = 0.364, as given 
by CELs, holds for ORLs too, then this brings the total ORL 
abundance up to O/H = 12.53 x 10"", i.e. 33 per cent higher. We 
adopt this latter value. 

NGC 3918: the C ii A4267, C ill AA4187, 4650 lines are detected; 

- 1 44 X 10"" is derived, as an intensity weighted mean 
from the AA4187, 4650 lines. Using the [Fe in] A4702 line, we 
estimate that 67 per cent of the detected A4658 intensity is due to 
C"+, i.e. /(C IV A4658) = 0.1093; hence, C"+/H+ = 2.62x 10"^ 
An ORL abundance is not available for C^, so it was estimated 
usms the Paper 1 CEL ratio of C+/C^+ = 0.155 ( given the small 
ionic concentration of C^, 9 per cent maximum, the error intro- 
duced should be negligible); adding this to the remaining three 
ionic stages we find, C/H = 7.50 x 10"" (if instead we correct for 
C""" using an ic/(C) = 1.09, we arrive at a value which is smaller by 
only 3 per cent). 

From an intensity weighted mean of four N^^/H^ ratios (Ta- 
ble 4), N^+/H+ = 1.42 X 10"" is found; from the N ill A4379 line, 
N^+/H+ = 6.41 X 10"'; all of the measured intensity of the A4606 
line is attributed to N IV (no contribution from N II A4607.2 is 
expected, since the stronger predicted N ii A4601 line is absent); 
therefore, N"+/H+ = 1.90 x lO"'^. We account for the missms N+ 
using the Paper 1 CEL ratio of N^/N^^ = 0.203; summing all four 
ionic stages, we have N/H = 2.56 x 10"". For this PN, we find that 
just as is the case for relative ionic fractions of carbon, the CEL 
N^+/N^+ ratio (= 1.22) is not equal to the ORL N^+/N^+ ratio (= 
0.45). 

We have derived 0^^/H^ and 0"^/H^ abundances from 
five O II multiplets (Table 6) and the O IV A4632 line (Table 7), 
respectively; 0^+/H+ = 5.36 x 10"" and 0"+/H+ = 4.74 x 10"^ 
We account for unseen 0+ and 0^+, using the Paper 1 CEL ratios 
of 0+/0^+ = 0.086 and 0^+/0^+ = 1.21; hence, summing all four 
ionic stages we have O/H = 12.32 x 10"". 

NGC 2022: the C ii A4267, C iii AA4069, 4187 and A4650 Unes as 
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well as the C IV A4658 line are detected; C-^+/H+ = 3.10 x 10"'' 
is the intensity weighted mean from the AA4187, 4650 lines only. 
We attribute all of the A4658 intensity to C iv, since other [Fe ill] 
lines are not present; furthermore this is a very high excitation ob- 
ject (He'^"'"/He = 0.88) and we expect a substantial concentration 
of €"+; we find C*+/H+ = 2.62 x 10"*. Using icf(C) = 1.02, to 
correct for unseen C"*" and summing the three higher ionic stages, 
we obtain C/H = 14.84 x 10"^. 

We derive N^+ZH"*" = 4.38 x 10"'' from an intensity-weighted 
mean of three N2+/H+ ratios (Table 4); from N III A4379, N''*+/H+ 
= 2.08 X 10"" and from N IV A4606, N''+/H+ = 3.68 x 10"^ (Ta- 
ble 5; as with NGC 3918, no contribution by N II A4607.2 is fotmd 
for the A4606 line). We correct for unseen N"*" using the Paper I 
CEL ratio N+/N'^+ = 0.013 and sum a total of four ionic stages to 
obtain N/H = 6.88 x 10"". 

From four O II multiplets (Table 6) and the O IV A4632 line 
(Table 7), respectively, we derive 0^+/H+ = 13.01 x 10"" and 
0*+/R+ = 1.10 X 10""; we also adopt 0^+/H+ = 3.62 x lO"", 
as derived by LD93a from the O III V8 multiplet at 3265 A. 
The unseen is accounted for using the Paper I CEL ratio of 
0+/0^+ = 0.020; hence we obtain 0/H = 17.99 x lO"". 

NGC 6818: the C ii A4267, C III AA4069, 4187 and A4650 
and C IV A4658 lines are detected; C^+/H+ = 1.37 x lO"", is 
the intensity-weighted mean from the AA4187, 4650 lines only 
(Table 2). Using exactly the same argument as in the case of 
NGC 2022, we attribute the A4658 line to C*+ only; we obtain 
C-*+/H+ = 4.47 X lO"--'. If we assume that the ratio C+/C^+ = 
0.129 derived from the CELs in Paper I is also valid for the ORLs, 
we deduce a total abvmdance of C/H = 7.71 x 10~*. 

The N II AA5676, 5679 lines from multiplet V3 are detected; 
they yield N^+/H+ = 4.27 x 10"" (Table 4); from N III A4379 we 
find N'^+/H+ = 8.87 x 10"^, while from N IV A4606 we deduce 
N''+/H+ = 1.46 X 10"^ (Table 5); the unseen N+ is corrected for 
using the Paper I CEL ratio N+/N^+ = 0.241. Summing a total of 
four ionic stages, N/H = 6.33 x 10~" is obtained. 

We have derived 0^+/H+ = 7.33 x 10"" and 0"+/H+ = 
6.41 X lO"'^ from four O II multiplets (Table 6) and the O IV 
A4632 line (Table 7), respectively; we further adopt O^'^/K^ 
= 2.37 X 10"", as derived by LD93a from the O ill V8 line at 
3265 A. The unseen O"'" is accounted for using the Paper I CEL 
ratio of 0+/0^+ = 0.090; hence, 0/H = 11.00 x 10"". 

NGC 3132: only C ii A4267 is detected. If we then use the stan- 
dard KB94 icf to account for unseen C^+, icf(C) = 0/0^+ = 1.93, 
and calculate the carbon abundance as C/H = icf(C) x C'^^/H^, we 
then find C/H = 12.74 x I0~". If instead, we assvune that the ratio 
q3+/q2+ ^ 0. 1 1 1 given by the CELs in Paper I is also valid for the 
ORLs and use it to estimate the ORL C'^'^flV' abundance, we find, 
C/H = 13.86 X 10"", i.e. 9 per cent larger. However, from observa- 
tions of 7 PNe it is found that, in general, ORL C*+/C^+ CEL 
C^+/C2+; hence, we adopt the former C/H value. 

From three N ii ORL lines and the N ill A4379 line we ob- 
tain, N2+/H+ = 3.40 X 10"" and N^+/H+ = 2.47 x 10"^ respec- 
tively. We correct for the unseen N^ using the Paper I CEL ratio 
N+/N^+ = 1.296 and sum a total of three ionic stages to find N/H 
= 8.05 X 10"". 

From the O II V 1 and V 10 multiplets (Table 6) we have 
derived O^'^/H^ = 8.15 x 10""; we account for unseen O'^^ using 
an icf{0) = 1.02 from Eq. A2; we account for unseen 0+ using the 
Paper I CEL ratio of 0+/0^+ = 0.894, hence, 0/H = 15.74 x 10"". 
The concentration of O"*" in this PN is about 46 per cent and it may 



be that some uncertainty has been introduced by assuming that the 
CEL 0+/0^+ ratio equals the ORL 0+/0^+ ratio. 

NGC 2440: the C ii A4267, C iii AA4187, 4650 and C iv A4658 
lines are detected; we attribute the A4658 intensity solely to C"+; 
thus C^+/H+ = 1.66 X 10"" and C"+/H+ = 1.82 x 10"" (Table 2). 
Using the KB94 icf{C)=l + 0+/0^+ = 1.28 to account for un- 
seen C+ and summing the three higher ionic stages, we find C/H = 
10.21 X 10"". 

Four N II lines are detected (Table 4); their intensity weighted 
mean yields N^+/H+ = 4.77 x 10""; from N ill A4379 we find 
N*+/H+ = 1.43 X 10""; both N iv AA4606, 4707 lines are de- 
tected, yielding a mean N"+/H+ = 1.70 X 10"". NoN+ abundance 
is available either from ORLs or CELs (our wavelength coverage 
did not include the [N II] lines). Thus, we adopt the KB94 CEL 
ratio of N+ZN^"*" = 0.475, and use it to correct for the vmseen N"*"; 
summing four ionic stages we obtain N/H = 10.17 x 10"". 

From four O II multiplets and the O IV A4632 line, re- 
spectively, we derive 0^+/li+ = 5.23 x 10"" and 0"+/H+ = 
1.01 X 10""; we further adopt 0-^+/H+ = 2.56 x 10"", as derived 
by Liu & Danziger (1993) from the O III V8 line at 3265 A, 
together with the CEL ratio 0+/0^+ = 0.285 from the same 
authors (as reported for their PA = 270" slit); we assume that this 
ratio is valid for ORLs too and hence derive O/H = 10.29 x 10"". 

NGC 6302: only C II A4267 is reliably detected; using the same 
icf scheme as for NGC 3132, yields icf (C) = 2.72, from Eq. Al we 
obtain C/H = 5.45 x 10"". 

Three N II lines are detected (Table 4); their intensity weighted 
mean yields N^+/H+ = 5.64 x 10""; N ill A4379 yields N*+/H+ = 
8.07 X 10""; N IV A4606 yields N"+/H+ = 2.28 x 10"" (no con- 
tribution is expected from N II A4607). N"'" is accounted for using 
the CEL ratio N"'"/N^"'" = 0.435; sunmiing four ionic stages we find 
N/H = 11.18 X 10"". 

We obtain 0^+/H+ = 3.28 x 10"" from the O ii V 1 multiplet; 
0"+/H+ = 7.84 X 10""^ is derived from O IV A4632; we assume 
that the Paper I CEL ratios of 0+/0^+ = 0.092 and 0^+/0^+ 
= 0.895 are valid for ORLs as well, and deduce 0/H = 7.30 x 10"". 

IC4406: C ii A4267, C iii AA4187, 4650 are detected (Table 2); 

the He^^ concentration is only 6 per cent, so it is assumed that the 
Q'i+/l{+ abundance is negligible. Using Eq. Al, with icf(C) =1.31 
to account for unseen C+, we find C/H = 1 1.27 x 10"" (Table 2). 

From the sole detected N ii A4630 line, N^+/H+ = 
1.90 X 10""; from N ill A4379, N^+/H+ = 7.31 x 10"^ N+ is ac- 
counted for using the CEL ratio N+/N^+ = 0.536; summing three 
ionic stages we find N/H = 3.65 x 10"" (Table 5). 

We obtain 0^+/H+ = 7.06 x 10"" and use an icf(0) = 1.04 to 
correct for unseen 0^+; the missing O^ is corrected for using the 
Paper I CEL ratio of O'^/O^^ = 0.471 (some error may be intro- 
duced since in this object, the O""" concentration is about 3 1 per cent 
of all oxygen). Hence, from Eq. A2 we obtain O/H = 10.80 x 10"" . 

1C4191: C II A4267, C III AA4187, 4650 are detected (Table 2); 
C'^+/H+ = 1.07 X 10"" is found. No significant C"+ is present. 
Using Eq. Al with icf{C) = 1.04 to account for unseen C"*", we find 
C/H = 6.48 X 10"". 

Six N II lines have been detected (Table 4); an intensity 
weighted mean yields N^+/H+ = 5.00 x 10""; from N III A4379, 
we derive N^+/H+ = 1. 15 x 10"" (Table 5). N+ is accounted for 
using the Paper I CEL ratio of N+/N^+ = 0.393, where N^+/H+ 
= 2.45 X 10"^ from the [N III] 57-/im line, as given by Liu et al. 
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(2001) and N+/H+ from our optical observations; hence, summing 
three ionic stages we obtain N/H = 8.12 x lO""*. 

We derive 0^+/H+ = 12.92 x 10"* from four O ii multiplets 
(Table 6) and use an icf(0) = 1 .07 to correct for unseen O^"*" ; the 
missing is corrected for using the CEL ratio of 0+/0'^+ = 
0.041 (the error introduced in this case should be negligible, since 
the O"*" concentration in this PN is only about 4 per cent). Hence 
we deduce 0/H = 14.39 x 10"*. 

My Cn 18: only C II A4267 is seen; no significant amounts of C''+ 
or C''"'" are expected; using Eq. Al with ic/(C) = 1.59 to account 
for unseen C+, we find C/H = 6.04 x 10"*. 

Three N II lines of multiplet V5 have been detected (Table 4); 
from an intensity weighted mean N^^/H^ = 20.35 x 10^* is de- 
rived. Apart from optical CELs no other information on N exists 
for this PN and the standard KB 94 icf scheme cannot account 
for unseen ionization stages. Therefore, as a lower limit we adopt 
N/H > 20.35 X 10"*. 

We derive 0^+/H+ = 6.43 x 10"* from the O II V 1 AA4649, 
4650 lines emitted by this low excitation nebula; unseen is 
corrected for using the Paper I CEL ratio of 0+/0^+ = 0.59 
(0+ represents a significant amount of the total oxygen, about 
37 per cent). The resulting 0/H = 10.24 x 10"* may be somewhat 
uncertain. 

SMCN87: C ii A4267, C iii AA4187, 4650 are detected (Table 3); 
C*+/H+ = 6.75 X 10"^ is found. No C*+ is present; using Eq. Al 
with Jc/(C) = 1.02 to account for the unseen only, we find C/H 
= 7.55 X 10"*. If instead we assume that the ratio of C+/C^+ = 
0.095 derived from CELs in Paper I is also valid for the ORLs, we 
arrive at the same total C/H value. 

We derive 0^+/H+ = 2.95 x 10"* from the O II V 1 multiplet 
(Table 3); unseen O"'" is corrected for using the Paper 1 CEL 
ratio 0+/0^+ = 0.018 (this is a medivun excitation PN and the 
concentration of 0+ is only about 2 per cent); we obtain 0/H = 
3.00 X 10"*. 

LMC N66: this is a highly ionized nebula reminiscent of the galac- 
tic Type-I PN NGC 6302; no ORLs of carbon are detected, con- 
sistent with the carbon-poor, nitrogen-rich nature deduced for this 
nebula from the CEL analysis of Paper 1, hence no ORL C/H value 
is available for this object. 

From an intensity weighted mean of O ii AA4649, 4650 
(Table 3), we derive 0^+/li+ = 8.53 x 10"*; we adopt ic/(0) = 
3.77 [from eq. (A7) of KB94 coupled with our Paper I CEL ionic 
nitrogen abundances for this PN]; we account for unseen O"'" using 
the CEL ratio 0+/0^+ = 0.076 and obtain 0/H = 3.46 x 10"^. 

LMC N141: C II A4267, C ill AA4187, 4650 are detected (Table 3); 
C^+/H+ = 5.82 x 10"^ is found. No C*+ is present; using Eq. Al, 
with !c/(C) = 1.02 to account for the unseen C+, we find C/H = 
7.54 X 10"*. If instead, we assume that the ratio C+/C^+ = 0.203 
derived from the CELs in Paper I is also valid for the ORLs, we 
obtain C/H = 8.75 x 10"*. In view of the error this may introduce, 
we adopt the former value. 

We derive 0^+/H+ = 4.96 x 10"* from the O II V 1 and V 10 
multiplets; imseen O""" is corrected for using the Paper I CEL ratio 
of 0+/0^+ = 0.023; we obtain O/H = 5.08 x 10"*. 

Nitrogen ORLs have not been detected from either SMCN87, 
LMCN66 or LMCN141, hence no inference can be made about 
their N abundances from optical recombination lines. 



APPENDIX B: A COMPARISON OF OBSERVED AND 
PREDICTED INTENSITIES OF O ii ORLS 

In Table B13 we present a comparison of the dereddened O II line 
intensities detected from all nebulae, against the predicted intensi- 
ties from recombination theory. The comparison is relative to the 
strongest expected line within each multiplet for the 3-3 transi- 
tions, while for the collection of 3d-4f transitions it is relative to the 
strongest expected 3d-4f line at 4089.3 A. The bracketed figures 
are the estimated absolute errors in the values, arising from the line 
profile fitting method only, and do not include any possible system- 
atic errors (for instance, those arising from the calibration process 
or the corrections for atmospheric and interstellar extinction). In the 
case of some of the 4f-3d Unes, additional, weaker, lines from the 
same multiplet contribute, as Usted in Table 4(a) of LSBC. The /obs 
values listed for some of the 4f-3d lines in Table B13 therefore in- 
clude contributions from weaker components of the multiplet, but 
the contributions from such blending lines were corrected for us- 
ing the theoretical branching ratios listed in Table 4(a) of LSBC in 
order to arrive at the /obs/^pred values in the final column of Ta- 
ble B13, which refer to the Usted multiplet component alone. 



Table B13. Comparison of observed and predicted relative intensi- 
ties of O II lines. 
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4129.32 


V19 


4-2 


1.70[.24] 


21. [3.0] 


4132.80 


V19 


2-4 


0.69[.15] 


1.0[0.2] 


4153.30 


V19 


4-6 


1.00 


1.0 


4156.53 


V19 


6-4 


1.23 [.20] 


7.7[1.3] 


4169.22 


V19 


6-6 


0.41[.07] 


1.2[0.2] 


3d-^f 










4083.90 


V48b 


6-8 


0.19[.03] 


0.7[0.1] 


4087.15 


V48c 


4-6 


0.18[.03] 


0.7[0.1] 


4089.29 


V48a 


10-12 


1.00 


1.0 


4275.55 


V67a 


8-10 


0.28[.05] 


0.4[0.1] 


4276.75 


V67b 


6-8 


0.25[.04] 


1.0[0.1] 


4277.43 


V67c 


2-4 


0.1 5 [.02] 


0.7 [0.1] 


4282.96 


V67c 


4-6 


0.08[.01] 


0.5[0.1] 


4283.73 


V67c 


4-4 


0.09[.01] 


0.9[0.1] 


4285.69 


V78b 


6-8 


0.11[.01] 


0.6[0.1] 


4288.82 


V53c 


2-4 


0.05[.02] 


0.5[0.2] 


4291.25 


V55 


6-8 


0.15[.04] 


0.9[0.2] 


4292.21 


V78c 


6-6 


0.08[.03] 


0.9[0.3] 


4294.78 


V53b 


4-6 


0.28[.03] 


1.0[0.1] 


4303.83 


V53a 


6-8 


0.52[.09] 


1.1 [0.2] 


4307.23 


V53b 


2^ 


0.08[.02] 


0.7[0.2] 


4313.44 


V78a 


8-10 


0.08[.02] 


0.7[0.2] 


4353.59 


V76c 


6-8 


0.09[.02] 


0.9[0.2] 


4357.25 


V63a 


6-8 


0.16[.02] 


2.7[0.3] 


4466.42 


V86b 


4-6 


0.23[.06] 


2.3[0.6] 


4477.90 


V88 


4-6 


0.13[.03] 


1.4[0.3] 


4489.49 


V86b 


2-4 


0.05[.02] 


0.7[0.3] 


4491.23 


V86a 


4-6 


0.21 [.02] 


1.5[0.1] 


4609.44 


V92a 


6-8 


0.20[.23] 


0.5 [0.4] 


4669.27 


V89b 


4-6 


0.11 [.02] 


2.8[0.5] 



A deep survey of heavy element lines in PNe - 

Table B13. -continued 



Ao(A) Mult. gl-gn /obs /obs//pred 



NGC5882 

3s-3p 



4638.86 


VI 


2-4 


0.47[.051 


2.2[0.11 


4641.81 


VI 


4-6 


0.70[.05] 


1.3[0.1] 


4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


0.33[.02] 


1.6[0.1] 


4661.63 


VI 


4-4 


0.33[.02] 


1.2[0.1] 


4673.73 


VI 


4-2 


0.05[.01] 


1.2[0.3] 


4676.24 


VI 


6-6 


0.25[.02] 


1.1[0.1] 


4317.14 


V2 


2-4 


0.70[.20] 


1.8[0.5] 


4319.63 


V2 


4-6 


0.37[.13] 


0.9[0.3] 


4345.56 


V2 


4-2 


1.21 [.50] 


3.0[1.3] 


4349.43 


V2 


6-6 


1.00 


1.0 


4414.90 


V5 


4-6 


1 .00 


1.0 


4416.97 


V5 


2-4 


0.91 [.22] 


1.6[0.4] 


op— 3a 










AflfiQ SO 


VI n 

V lU 


H — D 


U.'f7[. lOJ 


n 7rft '71 


4079 ^f\ 


V 1\J 


D — 


ft RTF 1 A1 


1 '^rft 91 


H-U / J . Ot) 


VIO 


o — iU 


1 00 


1 


4085 1 1 


VIO 


o — o 


n 1 1 r ClA^ 


Qrn 


4092 93 


VIO 


8-8 


osr 071 


9[0 2] 


4121.46 


V19 


2-2 


0.30[.23] 


0.8[0.6] 


4132.80 


V19 


2-4 


0.56[.08] 


0.8[0.1] 


4153.30 


V19 


4-6 


1.00 


1.0 


4156.53 


V19 


6-4 


0.76[.05] 


4.8[0.3] 


4169.22 


V19 


6-6 


0.29[.06] 


0.9[0.2] 


4110.78 


V20 


4-2 


0.39[.07] 


14[0.3] 


4119.22 


V20 


6-8 


1.00 


1.0 


4890.86 


V28 


4-2 


0.13[.07] 


0.5[0.3] 


4906.83 


V28 


4-4 


0.55[.17] 


0.9[0.3] 


4924.53 


V28 


4-6 


1.00 


1.0 


3d^f 










4083.90 


V48b 


6-8 


0.15[.05] 


0.5 [0.2] 


4087.15 


V48c 


4-6 


0.25[.04] 


0.9[0.1] 


4089.29 


V48a 


10-12 


1.00 


1.0 


4275.55 


V67a 


8-10 


0.40[.07] 


0.6[0.1] 


4276.75 


V67b 


6-8 


0.16[.07] 


0.7[0.3] 


4277.43 


V67c 


2-4 


0.18[.04] 


0.8[0.2] 


4291.25 


V55 


6-8 


0.1 5 [.03] 


0.9[0.2] 


4292.21 


V18c 


6-6 


0.05[.03] 


0.6[0.4] 


4294.78 


V53b 


4-6 


0.11 [.02] 


0.4[0.1] 


4303.83 


V53a 


6-8 


0.39[.04] 


0.8[0.1] 


4353.59 


V76c 


6-8 


0.10[.01] 


1.0[0.1] 


4366.53 


V75a 


6-8 


0.42[.05] 


0.6[0.1] 


4466.42 


V86b 


4-6 


0.21[.05] 


2.3[0.6] 


4491.23 


V86a 


4-6 


0.14[.03] 


1.0[0.2] 


4609.44 


V92a 


6-8 


0.42[.08] 


1.0[0.2] 
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Table B13. -continued 



Table B13. -continued 





Mult. 




-^obs 


T t IT 

-*obs'-*pred 




Mult. 


gl-gu 


-^obs 


T t IT 

-'obs'-'pred 






NGC5315 








NGC3918 




3s-3p 










3d-^f 










4638.86 


VI 


2-4 


0.49[.04] 


2.3[0.2] 


4083.90 


V48b 


6-8 


0.50[.18] 


1.7[0.6] 


4641.81 


VI 


4-6 


0.66[.05] 


1.3[0.1] 


4087.15 


V48c 


4-6 


0.43[.17] 


1.6[0.6] 


4649.13 


VI 


6-8 


1.00 


1.0 


4089.29 


V48a 


10-12 


1.00 


1.0 


4650.84 


VI 


2-2 


0.25[.03] 


1.2[0.1] 


4275.55 


V67a 


8-10 


1.15 


0.9 


4661.63 


VI 


4^ 


0.25[.02] 


0.9[0.1] 


4291.25 


V55 


6-8 


0.44 


1.7 


4673.73 


VI 


4-2 


0.04[.01] 


1.0[0.3] 


4294.78 


V53b 


4-6 


0.45 


1.6 


4676.24 


VI 


6-6 


0.22[.01] 


1.0[0.1] 


4303.83 


V53a 


6-8 


0.50 


1.1 


4317.14 


V2 


l-A 


1.00 


1.0 


4466.42 


V86b 


4-6 


0.43 


4.3 


4319.63 


V2 


4-6 


0.71[.10] 


0.7[0.1] 


4610.20 


V92c 


4-6 


0.67 


1.2 


4325.76 


V2 


2-2 


0.28[.04] 


1.5[0.2] 


4669.27 


V89b 


4-6 


0.22 


5.5 


4345.56 


V2 


4-2 


0.92[.20] 


1.0[0.3] 






NGC 2022 




4349.43 


V2 


6-6 


1.07[.20] 


0.5[0.1] 


3s-3p 










4366.89 


V2 


6-4 


0.85[.40] 


0.8[0.3] 


4638.86 


VI 


2-A 


0.32[.ll] 


1.5[0.5] 


4414.90 


V5 


4-6 


1.00 


1.0 


4641.81 


VI 


4-6 


1.10[.16] 


2.1 [0.3] 


4416.97 


V5 


2-4 


0.86[.50] 


1.5[0.9] 


4649.13 


VI 


6-8 


1.00 


1.0 


4452.37 


V5 


4^ 


0.20[.ll] 


1.8[1.0] 


4650.84 


VI 


2-2 


0.15[.21] 


0.7[1.0] 


3p-3d 










4661.63 


VI 


4-4 


0.33[.06] 


1.2[0.2] 


4069.62 


VIO 


4-6 


3.02[2.2] 


4.1 [3.0] 


4676.24 


VI 


6-6 


0.23[.08] 


1.0[0.4] 


4072.16 


VIO 


6-8 


1.33[1.0] 


1.9[1.4] 


3p-3d 










4075.86 


VIO 


8-10 


1.00 


1.0 


4072.16 


VIO 


6-8 


1.07[.17] 


1.6[0.3] 


4085.11 


VIO 


6-6 


0.31[.21] 


2.4[1.6] 


4075.86 


VIO 


8-10 


1.00 


1.0 


3d-^f 










3d^f 










4083.90 


V48b 


6-8 


0.27[.05] 


0.9[0.2] 


4089.29 


V48a 


10-12 


1.00 


1.0 


4087.15 


V48c 


4-6 


0.35[.04] 


1.3[0.2] 


4276.75 


V67b 


6-8 


0.92 


1.0 


4089.29 


V48a 


10-12 


1.00 


1.0 


4477.90 


V88 


4-6 


0.54 


6.0 


4275.55 


V67a 


8-10 


1.17[.20] 


1.0[0.2] 


4609.44 


V92a 


6-8 


0.29 


0.7 


4303.82 


V53a 


6-8 


0.37[.06] 


0.8[0.1] 






IC4406 




4609.44 


V92a 


6-8 


0.53[.04] 


1.2[0.1] 


3s-3p 














NGC3918 




4638.86 


VI 


2-A 


0.64[.ll] 


3.1 [0.5] 


3s-3p 










4641.81 


VI 


4-6 


0.93[.12] 


1.8[0.2] 


4638.86 


VI 


2-4 


1.04[.23] 


5.0[1.1] 


4649.13 


VI 


6-8 


1.00 


1.0 


4641.81 


VI 


4-6 


1.03[.40] 


1.9[0.7] 


4650.84 


VI 


2-2 


0.46[.09] 


2.2[0.4] 


4649.13 


VI 


6-8 


1.00 


1.0 


4661.63 


VI 


4-4 


0.35[.09] 


1.3[0.3] 


4650.84 


VI 


2-2 


0.26[.07] 


1.3[0.4] 


4676.24 


VI 


6-6 


0.17[.10] 


0.8[0.5] 


4661.63 


VI 


4^ 


0.35[.04] 


1.3[0.2] 


3p-3d 










4673.73 


VI 


4-2 


0.15[.03] 


3.5[0.8] 


4069.62 


VIO 


2-A 


1.62[.60] 


2.2[0.8] 


4676.24 


VI 


6-6 


0.31 [.02] 


1.4[0.1] 


4072.16 


VIO 


6-8 


1.49[.51] 


2.2[0.8] 


4317.14 


V2 


2-i 


1.00 


1.0 


4075.86 


VIO 


8-10 


1.00 


1.0 


4319.63 


V2 


4-6 


041 


04 






NGC 6818 




4325.76 


V2 


2-2 


0.34 


1.9 


3s-3p 










4349.43 


V2 


6-6 


1.53 


0.7 


4638.86 


VI 


2-4 


1.19[.30] 


5.7[1.4] 


4414.90 


V5 


4-6 


1.00 


1.0 


4641.81 


VI 


4-6 


2.73[.41] 


5.2[0.8] 


4416.97 


V5 


2-4 


0.56 


1.0 


4649.13 


VI 


6-8 


1.00 


1.0 


4452.37 


V5 


4-4 


1.26 


11. 


4650.84 


VI 


2-2 


0.29[.ll] 


1.4[0.6] 


3p-3d 










4661.63 


VI 


4-4 


0.47[.07] 


1.8[0.3] 


4069.62 


VIO 


4-6 


8.11[4.1] 


11 [5.6] 


3p-3d 










4072.16 


VIO 


6-8 


3.62[1.8] 


5.3[2.6] 


4069.78 


VIO 


2-4 


0.30[.26] 


0.4[0.4] 


4075.86 


VIO 


8-10 


1.00 


1.0 


4072.16 


VIO 


6-8 


0.40[.14] 


0.6[0.2] 


4085.11 


VIO 


6-6 


0.30[.22] 


2.3[1.7] 


4075.86 


VIO 


8-10 


1.00 


1.0 


4121.46 


V19 


2-2 


3.54 


9.8 


4085.11 


VIO 


6-6 


0.16[.06] 


1.2[0.5] 


4132.80 


V19 


2-A 


0.69 


1.0 






NGC 2440 




4153.30 


V19 


4-6 


1.00 


1.0 


3s-3p 










4156.53 


V19 


e-A 


1.42 


8.9 


4638.86 


VI 


2-4 


2.85[.56] 


14.[2.6] 


4169.22 


V19 


6-6 


0.18 


0.5 


4641.81 


VI 


4-6 


6.54[.37] 


9.4[0.7] 












4649.13 


VI 


6-8 


1.00 


1.0 












4650.84 


VI 


2-2 


0.20[.26] 


1.0[1.2] 












4661.63 


VI 


4-4 


0.35[.06] 


1.3[0.2] 












4676.24 


VI 


6-6 


0.71 [.08] 


3.2[0.4] 












4414.90 


V5 


4-6 


1.00 


1.0 












4416.97 


V5 


2-4 


0.80[.22] 


1.4[0.4] 












4452.37 


V5 


4-4 


1.70[.37] 


15.[3.3] 



Table B13. -continued 



A deep survey of heavy element lines in PNe - 

Table B13. -continued 



Ao(A) Mult. gl-gu /obs -fobs/-''pred 



IC4191 

3s-3p 

entire nebula 



4638.86 


VI 


2-4 


0.33[.06] 


1.6[0.3] 


4641.81 


VI 


4-6 


0.56[.06] 


1.1[0.1] 


4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


0.27[.05] 


1.3[0.3] 


4661.63 


VI 


4^ 


0.27[.02] 


1.0[0.1] 


4676.24 


VI 


6-6 


0.22[.02] 


1.0[0.1] 






fixed sUt 




4638.86 


VI 


2-4 


0.34[.04] 


1.6[0.2] 


4641.81 


VI 


4-6 


0.28[.02] 


0.5[0.04] 


4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


0.24[.03] 


1.2[0.2] 


4661.63 


VI 


4^ 


0.20[.03] 


0.7[0.1] 


4676.24 


VI 


6-6 


0.1 5 [.02] 


0.7[0.1] 


4317.14 


V2 


2-4 


1.00 


1.0 


4319.63 


V2 


4-6 


0.82[.25] 


0.8[0.3] 


4345.56 


V2 


4-2 


1.79[.64] 


1.9[0.7] 


4349.43 


V2 


6-6 


3.44[1.0] 


1.5[0.4] 


4366.89 


V2 


6^ 


1.58L36] 


1.6[0.4] 


4414.90 


V5 


4-6 


1.00 


1.0 


4416.97 


V5 


2-4 


0.81 [.08] 


1.4[0.1] 






entire nebula 




4414.90 


V5 


4-6 


1.00 


1.0 


4416.97 


V5 


2-4 


0.70[.19] 


1.3[0.4] 


4452.37 


V5 


4-4 


0.53[.14] 


4.8[1.3] 


3p-3d 














entire nebula 




4069.62 


VIO 


4-6 


2.45[.66] 


3.3[0.9] 


4072.16 


VIO 


6-8 


1.27[.36] 


1.8[0.5] 


4075.86 


VIO 


8-10 


1.00 


1.0 


4085.11 


VIO 


6-6 


0.17[.09] 


1.3[0.7] 






fixed slit 




4069.62 


VIO 


4-6 


1.74[.38] 


2.4[0.5] 


4072.16 


VIO 


6-8 


1.03[.25] 


1.5[0.4] 


4075.86 


VIO 


8-10 


1.00 


1.0 


4085.11 


VIO 


6-6 


0.08[.04] 


0.6[0.3] 


3d-4f 














entire nebula 




4083.90 


V48b 


6-8 


0.28[.13] 


1.0[0.5] 


4087.15 


V48c 


4-6 


0.39[.13] 


1.3[0.4] 


4089.29 


V48a 


10-12 


1.00 


1.0 


4275.55 


V67a 


8-10 


1.42[.28] 


1.1 [0.2] 


4282.96 


V67c 


4-6 


0.40[.13] 


0.9[0.3] 


4303.82 


V53a 


6-8 


0.63[.16] 


1.3[0.3] 


4609.44 


V92a 


6-8 


0.58[.10] 


1.3[0.2] 






fixed slit 




4083.90 


V48b 


6-8 


0.34[.07] 


1.2[0.3] 


4087.15 


V48c 


4-6 


0.33[.07] 


1.1 [0.2] 


4089.29 


V48a 


10-12 


1.00 


1.0 


4275.55 


V67a 


8-10 


1.18[.18] 


1.0[0.2] 


4303.82 


V53a 


6-8 


0.33[.07] 


0.7[0.2] 


4466.42 


V86b 


4-6 


0.19[.03] 


1.9[0.3] 


4609.44 


V92a 


6-8 


0.29[.08] 


0.7[0.1] 



Ao(A) Mult. gl-gu /obs lobJIpred 



NGC3132 

3s-3p 



4638.86 


VI 


2-4 


0.53 [.471 


2.6[2.2] 


4641.81 


VI 


4-6 


0.88[.58] 


1.7[1.1] 


4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


0.66[.38] 


3.2[1.8] 


^00 1 .03 


V 1 


4-4 


0.60[.34] 




4676.24 


VI 


6-6 


0.39[.23] 


1 7ri 11 

1 . / L 1 . i J 


3p-3d 










4069 62 


VIO 


2-4 


1.14[.78] 


1.5[1.0] 


4072.16 


VIO 


6-8 


041[.34] 


0.6[0.5] 


4075.86 


VIO 


8-10 


1.00 


1.0 


4085.11 


VIO 


6-6 


0.12[.03] 


0.9[0.2] 






NGC6302 




3s-3p 










4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


0.97 


4.7 


4661.63 


VI 


4-4 


1.14 


5.1 






MyCnlS 




3s-3p 










4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


0.37 


1.8 






LMC N66 




3s-3p 










4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


1.00 


4.8 






LMC N141 




3s-3p 










4638.86 


VI 


2-4 


0.20[.10] 


1.0[0.5] 


4641.81 


VI 


4-6 


0.95 [.28] 


1.8[0.5] 


4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


0.66[41] 


3.1[1.9] 


4661.63 


VI 


4-4 


041 [.14] 


1.6[0.6] 


4676.24 


VI 


6-6 


0.37[.09] 


1.7L0.4] 


3p-3d 










4069.62 


VIO 


2-4 


0.74[.36] 


1.0[0.5] 


4072.16 


VIO 


6-8 


0.49[.18] 


0.7[0.3] 


4075.86 


VIO 


8-10 


1.00 


1.0 






SMC N87 




3s-3p 










4638.86 


VI 


2-4 


1.07[.17] 


5.1L0.8] 


4641.81 


VI 


4-6 


1.06[.26] 


2.0[0.4] 


4649.13 


VI 


6-8 


1.00 


1.0 


4650.84 


VI 


2-2 


1.00[.15] 


4.8[0.7] 


4661.63 


VI 


4-4 


0.31 [.04] 


1.2[0.2] 


4676.24 


VI 


6-6 


0.32[.ll] 


1.4[0.5] 



